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ABSTRACT 


Studies  of  superconductors  as  a  function  of 
pressure  have.,  until  now,  been  concerned  mainly  with 
the  bulk  properties  of  the  metals.  A  study  of  the  shift 
of  some  of  the  microscopic  properties  of  strong-coupling 
Pb  with  pressure  has  been  undertaken  and  this  thesis 
reports  the  results  of  this  pressure  study.  The  method 
of  investigation  used  was  the  electron  tunneling  technique, 
where  the  pressure  production  was  achieved  by  suspending 
the  samples  in  a  solid  helium  bath.  The  energy  gap  and 
phonon  spectrum  of  Pb  as  reflected  in  the  junction  con¬ 
ductance  curves  were  studied  as  a  function  of  pressure 
for  the  first  time. 

The  average  value  of  dT  /dP  from  many  different 
junctions  is  in  agreement  with  published  results  for 
bulk  Pb .  Shifts  in  the  predominant  transverse  and 
longitudinal  phonon  groups  were  observed.  These  shifts 
occur  in  the  dispersive  region  of  the  two  main  phonon 
branches.  Gruneisen  q's  of  7(0^)  =  2052  and  y(co^)  =  3.32 
were  obtained  for  the  phonon  peak-  regions.,  and  these  can 
be  favourable  compared  with  results  from  thermal  expan¬ 
sion  work. 

The  "recombination  effect"  was  studied  and  found  to 
be  smaller  than  predicted,  and  perhaps  absent  altogether. 
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The  reduced  relation  ( Aq( T)/kTc)/( Aq( T)/kT  was 

studied  as  a  function  of  reduced  temperature,  and  pressure. 
The  results  indicate  that  2AQ(0)/kTc  is  pressure  dependent., 
in  contrast  to  accepted  speculations.  Also,  the  reduced 
energy  gap  A  (T)/A  (0)  as  a  function  of  reduced  temperature 
t,  follows  the  BCS  relation  to  within  a  few  percent;  even 
under  pressure. 

The  frequency  dependence  of  A,  as  defined  by  McMillan 

o 

(1968),  was  found  to  be  cc  l/<a)  >  in  agreement  with  theory. 

The  electron-phonon  coupling  constant  (a  (cd))  was 
found  to  decrease  under  increased  pressure.  This  reduction 
in  a  (cd)  results  in  a  shift  in  the  properties  of  strong¬ 
coupling  superconducting  Pb  toward  the  predictions  of  the 
BCS  weak-coupling  theory. 
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CHAPTER  I 

INTRODUCTION 

In  1911  Onnes  discovered  superconductivity  while 
measuring  the  electrical  resistance  of  mercury  wires  at 
low  temperature.  Since  then  the  state  has  been  found  in 
many  other  metals  and  alloys,  and  is  reversible  in  the 
sense  of  a  true  thermodynamic  phase  state. 

Experiments  by  Meissner  and  Ochsenfeld  (1933) 
established  that  a  superconductor  expels  all  magnetic 
flux  from  its  interior.  (We  limit  ourselves  to  a  dis¬ 
cussion  of  Type  1  superconductors.)  It  was  found  that 
by  applying  a  magnetic  field  greater  than  some  critical 
value,  H  ,  the  normal  resistive  state  could  be  restored. 
The  critical  field  was  temperature  dependent  below  the 
superconducting  transition  temperature,  T  ,  the  depend¬ 
ence  being  approximately  of  the  form 

Hc  =  V1  -  (T/T0)2). 

An  abrupt  jump  in  the  specific  heat  at  the  transition 
temperature  in  zero  magnetic  field,  also  indicated  that 
the  transition  was  a  phase  change  of  second  order. 

Phenomenological  theories  were  proposed  first  by 
Gorter  and  Casimir  (193^0  and  later  by  F.  and  H.  London 
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(193^-j  1935).  These  were  based  on  a  two-fluid,  model 
concept o  One  fluid,  the  normal  component,  consisted  of 
ordinary  electrons,  while  the  second  fluid  consisted  of 
the  super-electrons.  In  the  London  theory,  the  normal 
fluid  was  treated  using  classical  electrodynamics,  while 
the  concept  of  infinite  conductivity  and  perfect  dia¬ 
magnetism  were  incorporated  in  dealing  with  the  super¬ 
fluid  component. 

I . 1  Experiments  Leading  up  to  a  Microscopic  Theory 

In  1950  Frohlich  theorectically  predicted  that  the 

critical  temperature  should  depend  on  the  isotopic  mass 

1. 

of  a  superconductor  according  to  the  relation  T^M2  =  const. 
This  was  experimentally  found  to  be  the  case  independently 
by  Maxwell  (1950)  and  by  Reynolds  et  al  (1951) »  This 
result  established  the  fact  that  the  lattice  must  be 
involved  in  the  mechanism  for  superconductivity. 

Between  1950  and  1957 ^  many  different  experiments 
were  suggesting  simultaneously  that  the  superconducting 
state  exhibits  an  energy  gap  2A  in  the  single  particle 
density  of  states  function  near  the  Fermi  surface.  These 
experiments  included  specific  heat  measurements,  infra¬ 
red  transmission  experiments,  electromagnetic  radiation 
absorption  measurements,  and  ultrasonic  attenuation 
experiments „ 


. 
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These  results,  coupled  with  the  earlier  known  result 
that  most  superconductors  obey  the  same  laws  of  corres¬ 
ponding  states  to  a  good  degree,  formed  the  starting 
point  for  a  microscopic  theory. 

Frohlich  (1950)  proposed  a  theory  of  superconductivity 
based  on  an  attractive  electron-electron  interaction 
mediated  by  virtual  phonon  exchange.  The  electrons  were 
not  necessarily  from  states  of  opposite  momenta,  a 
refinement  later  proposed  by  Cooper  (1956).  In  detail 
Frohlich  established  the  criterion  for  superconductivity, 
namely  that  if  the  attractive  electron-electron  inter¬ 
action  dominates  the  repulsive  Coulomb  interaction  between 
the  electrons,  a  metal  becomes  a  superconductor.  However, 
due  to  his  perturbation  approach,  which  could  not  be 
carried  out  with  sufficient  accuracy,  the  quantitative 
features  of  the  theory  failed. 

Cooper  (1956)  showed  that  a  pair  of  electrons  of 
opposite  momentum  and  spin,  interacting  via  a  momentum 
dependent  potential  above  a  filled  Fermi  sea  of  non¬ 
interacting  electrons,  would  form  a  bound  pair.  The 
potential  must  be  attractive  in  the  region  of  the  Fermi 
surface  although  it  could  be  arbitrarily  weak.  Finally, 
Bardeen,  Cooper,  and  Schrieffer  (1957)*  henceforth 
(BCS),  succeeded  in  producing  a  successful  theory  which 
involved  macroscopic  occupation  of  such  pair  states. 
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Simultaneously,,  Bogolyubov  (1958)  in  Russia  formulated 
an  independent  theory  based  on  a  somewhat  different 
mathematical  approach. 

More  recently.,  a  sophisticated  description  of  the 
behaviour  of  superconductors  has  become  possible  through 
the  use  of  field  theoretic  Green’s  functions  techniques. 
These  calculations  make  use  of  information  regarding 
electron-phonon  coupling  strengths  and  explicit  phonon 
distributions . 

Because  of  the  dependence  of  superconductivity  on 
the  lattice*  studies  of  superconductors  subjected  to 
pressure  have  been  carried  on  for  some  time.  The  lit¬ 
erature  has  been  reviewed  by  Brandt  and  Ginzburg  (1965), 
by  Levy  and  Olsen  (1965)*  by  Olsen  and  Rohrer  (i960), 
and  by  Swenson  (i960)  in  his  "High  Pressure  Physics" 
survey. 

The  greater  part  of  the  work  done  to  date  on  study¬ 
ing  the  influence  of  pressure  on  superconductors  has 
been  in  determining  the  sign  and  magnitude  of  c)T  /dP 
and  dHc/dP.  Characteristically,  most  superconductors 
decrease  in  both  T  and  H  with  increasing  pressure. 

C  C  o 

This  is  the  case  in  Pb .  Bismuth  is  interesting  in  that 
it  is  not  a  superconductor  unless  a  pressure  of  from 
20,000  -  40,000  atmospheres  is  applied  to  it.  Thus  it 
is  felt  that  only  one  phase  of  Bi,  the  one  not  normally 
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found  at  room  temperature,  will  superconduct  at  low 
temperature . 

An  increase  in  the  critical  temperature  is  observed 
in  zirconium,  titanium,  and  a  small  number  of  binary 
alloys.  Thallium  shows  an  increase  in  T  to  about  2000 
bars  and  then  it  begins  to  decrease. 

There  seems  to  be  very  little  experimental  infor¬ 
mation  on  the  shift  in  microscopic  parameters  of  super¬ 
conductors  under  pressure.  The  only  published  data  are 
from  experiments  by  Hodder  and  Briscoe  (1966)  and  by 
Zavaritskii  (1967) ^  both  groups  worked  with  Pb .  Hodder 
and  Briscoe  have  subjected  Pb-Pb  tunnel  junctions  to 
mechanical  bending.  They  observe  a  shift  in  the  gap- 
parameter  Aq  which  they  try  to  relate  to  a  volume  depend¬ 
ence.  However,  their  method  is  very  inaccurate  because 
of  shear  strains  set  up  in  the  specimen.  More  recently, 
Zavaritskii  has  investigated  the  strain  dependence  of 
the  phonon  spectrum  in  Pb .  He  deposits  the  Pb  strips  of 
a  tunnel  junction  at  low  temperature  and  in  this  way 
achieves  maximum  lattice  distortion.  By  annealling  the 
junctions  to  successively  higher  temperatures  and  measur¬ 
ing  the  accompanying  phonon  distributions,  he  obtains  some 
information  about  the  behaviour  of  a  superconductor  with 
varying  amounts  of  strain  in  the  lattice. 
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It  was  decided  that  experiments  on  tunnel  junctions 
subjected  to  hydrostatic  pressure,,  would  be  attempted. 
Provided  the  technique  was  feasible.,  this  would  be  the 
best  way  to  measure  the  pressure  dependence  of  the  energy 
gap  parameter  A  .  Tunnel  junctions  made  from  Al-Al-oxide- 
Pb  were  chosen  with  the  object  of  studying  superconducting 
Pb .  Pb  is  highly  compressible,,  has  a  large  Tc  and  gap 
value ,  and  exhibits  strong  phonon  structure  in  its 
density  of  states  which  could  serve  as  a  means  of  measur¬ 
ing  phonon  shifts  under  pressure. 

The  next  Chapters  briefly  review  existing  super¬ 
conductivity  theory  and  its  application  to  strong¬ 
coupling  Pb .  A  discussion  of  the  cryostat  and  electronic 
apparatus  is  given,,  followed  by  the  laboratory  procedure. 
Results  are  presented  and  discussed  in  Chapter  VI. 
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CHAPTER  II 


THEORY  OF  SUPERCONDUCTIVITY 


I I . 1  The  Electron-Phonon  Interaction 

The  electron-phonon  interaction  plays  an  essential 
role  in  determining  the  physical  properties  of  a  metal. 
In  calculating  effects  due  to  interaction  between  the 
electron  and  phonon  systems,  the  many  body  Hamiltonian 
for  the  system  is  first  set  down  in  terms  of  ’’bare” 
particles  consisting  of  single  electrons  occupying 
Bloch  states,  and  quantized  vibrations  of  the  ionic 
lattice.  Next,  coupling  between  these  particles  leads 
to  a  set  of  "dressed"  particles  in  terms  of  which, 
properties  of  the  physical  system  can  be  derived. 

The  nature  and  strength  of  the  electron-phonon  coupling 
has  important  consequences  for  the  properties  of  metals 

Included  in  the  zero  order  Hamiltonian  H  is  a 

o 

term  representing  the  interaction  of  the  electrons  with 
the  ions  in  their  equilibrium  positions.  The  differ¬ 
ence  between  this  potential  and  the  full  electron-ion 
potential  remains  as  a  perturbation  which  is  even¬ 
tually  related  to  an  electron-phonon  coupling  parameter 


In  the  case  of  weak  coupling  between  the  electron 


'  ■ 


8 

and  phonon  systems,,  the  electron  excitations  are  long 
lived  and  their  energies  well  defined.  This  seems  to 
be  the  case  whether  the  interaction  is  retarded  or 
not . 

On  the  other  hand.,  when  the  interaction  is  strong, 
the  retardation  of  the  electron-phonon  interaction 
leads  to  damping  mechanisms  and  hence  particle  life¬ 
time  effects  when  considering  the  excitations  of  the 
electron  distribution.  (These  lifetime  effects  are 
not  present  if  an  instantaneous  interaction  is  allowed.) 
The  form  of  the  phonon  spectrum  as  well  as  the  electron 
distribution  (Fermi  surface  shape)  play  significant 
roles  in  determining  the  predictions  of  a  specific 
theoretical  calculation. 

It  is  generally  assumed  that  the  electron-phonon 
interaction  is  essential  in  bringing  about  super¬ 
conductivity  in  most  (if  not  all)  superconductors. 

The  electron-phonon  interaction  can  lead  to  electron- 
electron  coupling  in  the  following  processes: 


■ 
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either  (a)  an  electron  in  state  k  emits  a  phonon 
and  is  scattered  to  k  -  K  while  electron  k'  absorbs 

it  and  is  scattered  to  k1  +  K,  or  (b)  a  phonon  with 

^  ^ 

wave  vector  -q  is  emitted  by  k' .  By  second  order 
perturbation  theory ,  the  matrix  elements  for  the 
transitions  are 


* 


K,k'  +  K  |  V  |  k,k'\  =  — — kUk'  +X 

ek  "  ^  ek-K  + 


and 


K,kl  +  K  |  V  |  k,k^>  = 


•X* 

M  M 

k,  k-K  k',k'  +  K 

€,  -  [  6,  -rr  Ilcd  ] 

k  k-K  qJ 


II-l 


II-2 


where  the  electron  energy  is  measured  relative  to  the 
Fermi  surface,  ttcu  is  the  energy  of  the  exchanged  phonon 
and  K  is  its  wave-vector  (K  =  q  for  normal  processes) . 

M  is  the  matrix  element  for  the  electron-phonon  inter¬ 
action  o  Combining  these  under  the  assumption  that  the 
electron-phonon  matrix  elements  M  and  M  are  the  same 
in  the  energy  range  of  interest  gives: 


(k  -  K,k»  +  K  |  V  |  k,k^> 


2  I  M 


k,  k-K 


R. 


^  €k  ek-K 


J2  -  (fta)  )2 


H-3 


Clearly  for  the  energy  region  |  |  < 

this  interaction  is  attractive.  In  the  normal-super¬ 
conducting  transition,  only  processes  within  an  energy 
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kTc  of  the  Fermi  surface  are  important,  where  Tq  is 
the  transition  temperature.  Since  generally  kTo«fia> 
for  <Dq  of  the  order  of  the  Debye  frequency, 

■  ,  2 

V(K)  ^  ^  1  k,k-K  1  . 


q 


II-4 


A  second  possible  interaction  between  electrons  is 
the  Coulomb  repulsion,  which  can  be  written  in  terms 
of  a  screening  parameter  A  as 


n  (K)  = 
Coulv  ’ 


47re‘ 


2  2 
Am  YT 


II-5 


The  total  interaction  will  be  therefore  given  by 

V(K)  =  ~2  '  Mk,k-K  I  +  Z|7re2  .  II- 6 

A2+  K2 

If  V(K)  is  negative,  an  attractive  force  is  felt  between 
the  electrons. 


II. 2  The  Bardeen,  Cooper,  Schrieffer  Theory 

In  1957^  BCS  formulated  the  first  successful 
theory  of  superconductivity.  Cooper  (1956)  had  found 
earlier,  that  when  two  electrons  interact  with  an 
arbitrarily  small  attractive  interaction  V  above  a 


'  \ 
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filled  sea  of  non-interacting  electrons,  binding  results 
between  the  pair.  Such -a  bound  pair  is  called  a  Cooper 
pair  o 

The  BCS  theory  considers  a  system  in  which  a  large 
number  of  electrons  interact  via  the  potential  V  in 
pairs.  The  matrix  elements  which  take  a  pair  of  elec¬ 
trons  from  any  two  initial  states  to  any  other  two 
states  satisfying  momentum  conservation,,  alternate  in 
sign  because  of  Fermi-Dirac  statistics.  If  these  con¬ 
figurations  occurred  in  the  ground  state  with  roughly 
equal  weight,  the  net  interaction  would  be  small 0  BCS 
showed  that  a  coherent  low  energy  state  can  be  produced 
by  choosing  a  subset  of  configurations  between  which  the 
matrix  elements  are  negative.  This  subset  consists  of 
states  composed  completely  of  paired  states,  all  having 

-A 

the  same  momentum  k  +  kl  =  C.  The  greatest  lowering  in 

v 

energy  occurs  when  the  greatest  amount  of  pairing  is 
produced.  For  ground  state  pairing,  this  occurs  when 
C  =  0,  as  can  be  seen  from  the  following  diagram  by 
Cooper . 


. 
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The  possible  values  of  k  and  k'  lie  in  the  narrow  shell 
straddling  k^ .  Also  those  values  satisfying  k  +  k'  =  G 
lie  in  the  two  shaded  regions .  This  shows  that  the 
maximum  amount  of  pairing  is  obtained  when  C  =  0. 

Further  study  shows  that  exchange  energy  consid¬ 
erations  require  that  pair  electrons  are  of  opposite 
spin,  since  the  interaction  is  reduced  for  parallel 
spins . 


Thus  BCS  consider  a  reduced  problem  in  which  only 
configurations  where  states  are  occupied  in  pairs  k  l 
and  -klr  are  important.  They  introduce  a  reduced 
Hamiltonian  which  includes  only  the  interactions  be¬ 
tween  zero  momentum  pairs,  and  this  presumably  accounts 
for  most  of  the  condensation  energy.  A  Hartree-type 
trial  wave  function  ^(u^v^ 

2 

satisfies  the  pairing  condition.  Here  v^  is  the 

th  2 

probability  that  the  k  pair  is  occupied  and  u^  is 

the  probability  that  it  is  unoccupied. 

The  energy  corresponding  to  the  reduced  Hamiltonian 
is  calculated  treating  v^  and  u^  as  variational  para¬ 
meters.  A  minimum  value  in  the  condensation  energy  W 

0 

is  obtained  when  these  parameters  are  given  by: 


)  is  introduced,  which 


1(1  -  ek/Ek) 


n-7 


and 


id  +  ck/Ek- 


II-8 
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The  energies  are  given  by  the  dispersion  relation 


+  A 


k 


H-9 


where  is  the  ordinary  Bloch  energy  measured  from 
the  Fermi  surface,  and  A^  is  defined  by 


Z 

k 


ik,uk,vk' 


II- 10 


Substituting  for  u^.,  and  v^,  produces  a  self-consistency 
condition  given  by. 


Z 

k 1 


V,  ,  .A,  . 
kk !  k  * 


2E. 


k 
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This  is  referred  to  as  the  gap  equation.  The  energies 
E^  turn  out  to  be  the  single  particle-like  excitation 
energies  of  the  system.  BCS  neglect  anisotropic  effects 
and  assume  that  V-^t  or  V(K)  can  be  replaced  by  a  constant 
average  matrix  element  such  that 

v  =  <V(K))  Av  for  I  €  I  < 


=  0 


elsewhere . 


Here  oo  is  an  average  phonon  frequency  chosen  equal  to 
Ql 

the  Debye  frequency  od^.  BCS  also  assume  that  the  density 
of  Bloch  states  is  constant  and  equal  to  the  value  at  the 
Fermi  level  (N(0))  in  the  range  |  e  |  <  . 
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The  assumption  that  V  =  constant  leads 
and 


to  A, 
k 


_ q 

sinh(l/N(0)V) 


constant,, 


11-12 


Aq  is  called  the  gap  parameter.  Under  these  same  assump¬ 
tions  the  condensation  energy  becomes^ 


-2N(0)(»m  )2 

w  -  w  =  w  =  — - q - 

s  n  c  exp[2/N(0)V]  -  1 


II-13 


o 

Empirically,  W  is  of  the  order  of  N(0)(kT  )  and  in 

c  c 

general  kT  is  much  less  than  frcu  .  According  to  11-13 
c  q 

this  will  occur  if  N(0)V  «  1,  which  is  termed  the  weak- 
coupling  limit.  Then  we  have  the  simplified  expressions: 


P  "2 

W  =  -2N(0)(fico  )  exp[ - 

°  q  N(  0)  V 


] 


II- 14 


-1 

A  =  2fiou  exp[ - 

0  q  N(0)V 


] 
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Excitations  above  the  ground  state  occur  when  one 
or  more  of  the  ground  state  pairs  is  broken  up.,  or  when 
pairs  are  excited.  It  was  shown  by  BCS  that  the  excita¬ 
tion  energy  for  creating  one  unpaired  quasi-particle  in 
state  k  is  given  by  E^  =  ( -f  AQ)2.  However,  such 
excitations  can  only  be  created  in  pairs;  the  minimum 


excitation  energy  is  therefore  2Aq. 


It  was  also  shown 
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that  excited  pairs  can  be  formally  treated  as  a  pair  of 
excited  quasi-particles.  The  excitations  behave  like 
independent  fermions.  We  have,  therefore,  a  one-to-one 
correspondence  between  excited  states  and  the  Bloch-type 
excitations  in  the  normal  metal. 

One  can  use  this  fact  to  calculate  the  density  of 
excited  states  in  a  superconductor.  We  have  from  this 
one-to-one  relationship: 

Vek>dek  =  Ns(EddEk  XI-16 

or 

VEk)  =  Nn(ek)^y-  X1-^ 

K 

Differentiating  the  dispersion  relation  gives: 


W  = 


Nn(£k>  S 


k 


£,  +  A  — 
k  de 


k 
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For  a  constant  gap  parameter,  the  normalized  super¬ 
conducting  density  of  states  becomes: 


n  (Eb)  .  W 

w 


Re 


(E2  -  A2)^ 


11-19 


and  n  is  sharply  peaked  at  energies  given  by  E  =  +  A. 
s 
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II . 3  The  BCS  Theory  at  Finite  Temperature 

The  behaviour  of  the  superconductor  for  T  >  0°K  is 
obtained  by  BCS  by  minimizing  its  free  energy.  They 
obtain  expressions  very  similar  to  the  ground-state 
calculations.  In  particular,  the  quasi-particle 
energies  are  given  by: 

Ek =  (ek  +  A0(T)2h  u-20 

where  A-  (T)  is  now  the  temperature-dependent  energy- 
gap  0  The  density  of  states,  equation  11-19,  is  there¬ 
fore  temperature-dependent.  A  (T)  is  found  to  go  to 
zero  at  some  finite  temperature,  T^,  and  this  is  recog¬ 
nized  as  the  transition  temperature.  In  the  weak-coupling 
limit,  N(0)V  «  1,  the  following  relations  hold: 

kT  =  l.l4fia)  exp (  -  1/N(0)V)  11-21 

2A  (0)  =  3.528  kT  .  11-22 

The  reduced  energy  gap  Aq(T)/A  (0)  is  a  universal 

function  of  reduced  temperature  t  =  T/T  .  In  particular, 

0 

near  T  ,  the  following  expressions  hold: 

0 

Aq(T)  =  3.069  kTc(l  -  T/Tc)2  11-23 

or 

Ao(T)/Aq(0)  =  1.74  (1  -  T/T,)t  11-24 

A  complete  tabulation  of  Aq  and  other  parameters  of  the 
BCS  superconductor  has  been  given  by  Muhlschlegel  (1959). 
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It  Is  well  known  from  a  thermodynamical  treatment , 

that  the  transition  at  T  =  T  is  of  second  order,  exhibit- 

c  J 

ing  a  finite  discontinuity  in  specific  heat.  The  BCS 
theory  predicts  for  this: 


C 


es 


-  C 


en 


C 


en 


3  d(A/kTc)2 
2i t2  dt 


which  gives  the  universal  number: 


11-25 


C 


es 


-  C 


en 


=  1.43 


C 


en 
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II o 4  Tunneling  and  the  Superconducting  Density  of  States 

A  detailed  one-to-one  correspondence  between  tunnel 
junction  conductance  and  the  superconducting  density  of 
states  has  been  established  by  Bardeen  (1961)  and  by  Cohen., 
Falicov,  and  Phillips  (1962).  A  simplified  version  of  this 
will  be  outlined  briefly  in  the  present  section. 

A  tunnel  junction  consists  of  a  metal-insulator- 
metal  sandwich  formed  by  first  depositing  a  metal  layer 
onto  a  suitable  backing.,  allowing  its  surface  to  oxidize 
forming  the  insulator.,  and  then  cross-depositing  a  second 
metal  layer .  Provided  the  barrier  oxide  layer  is  suffi¬ 
ciently  thin,  electrons  can  penetrate  through  it  by 
quantum  mechanical  tunneling. 


■ 
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A  simple  model  for  the  tunnel  current  exists  for 
a  normal- superconducting  (n-s)  junction  (Shapiro  et  al, 
1962),,  that  is,  a  junction  in  which  one  metal  is  in 

the  normal  state  and  the  other  is  in  the  superconducting 

2 

state o  Assuming  a  tunneling  probability  |M| for  barrier 
penetration  which  is  to  be  taken  independent  of  tempera¬ 
ture  and  energy  for  a  bias  of  up  to  ^  30  millivolts, 
the  transition  probability  per  unit  time  for  an  electron 
in  an  occupied  state  k-^  tunneling  to  an  unoccupied 
state  kg  on  the  other  side  of  the  barrier  is: 

P  =  2ir  |  M  |  2f\(l  -  ff)Nf  11-27 

is  the  density  of  final  states  and  fn.  and  f^  are 
Fermi  functions  for  the  initial  and  final  states.  If  a 
voltage  V  is  now  impressed  across  the  junction.,  the  net- 
tunneling  current  will  be  the  sum  of  the  two  one-way 
currents.  Thus,  (gathering  all  the  constant  factors  into 
one) 

00 

i(V)  =  const  e  /(N1(E)f (E)Np(E  +  eV) [ 1  -  f(E  +  eV) ] 

-00 

-  [N1(E)f(E  +  eV)Ng(E  +  eV) ( 1  -  f(E))]}dE  11-28 

where  E  is  the  energy  measured  from  the  Fermi  level. 
Collecting 

00 

i(V)  =  const  e  /N, (E)Ng(E  +  eV)[f(E)  -  f(E  +  eV)JdE. 

11-29 


- 
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Recalling  equation  11-19*  the  density  of  states  for  a 
superconductor ,  with  energy  measured  from  the  Fermi  level, 
is  given  by  the  form 


Nt(o>) 


and  for  go 
(choice 
forgetting 


i  (V) 
ns v  ’ 


=  N(0)  Re  f  \ 

\  p  p  jU 

(cd  -  A  J 2 

=  eV  <  A,  Nrp(o>)  =  0.  Integrating  11-29 
of  go  =  0)  to  go  =  eV  for  a  n-s  junction, 
factors  of  two,  for  T  =  0 

=  const  e  (0)^  (0)  [  (eV) 2  -  A"-]2. 


n-30 


from 

and 


11-31 


Therefore  the  junction  conductance,  assuming  the 
normal  density  of  states  is  constant,  is 

di(V)  eV 

-  =  const  e  N  '(O)Np(O) - ^ - 75— r  11-32 

dv  1  d  ((eV)d  -  A  ; 2 


and  the  normalized  conductance  is 


(di/dV)ng  Re  |  eV  |  NT(cu) 

<J  ==  —  p  p  3~  =  •  1 1  ”  3  3 

(di/dV)nn  ( (eV;  -  A")2  N(0) 

Thus  the  normalized  conductance  gives  a  direct  measure 
of  the  normalized  density  of  states  in  a  superconductor. 

Coherence  effects  associated  with  the  interdependence 

_x  t  -x 

of  an  electron  in  state  kf  with  an  electron  in  state  -k  l 
often  arise  when  considering  single  particle  processes  in 
a  superconductor.  These  effects  do  not  occur,  however. 
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when  the  one-particle  tunneling  process  between  a 
normal  and  superconducting  metal  is  considered.  Consider 
the  illustrations  in  figure  13  due  to  Schrieffer.  An 
electron  in  k  f  beneath  the  Fermi  surface  in  the  normal 
metal  tunnels  through  the  oxide  to  state  k’l'  above  the 
Fermi  surface  in  the  superconductor.  In  the  process  a 
hole  is  left  behind  giving  an  excitation  energy  e  =  |  e 

IX 

for  this  metal.  In  addition a  quasi-particle  Is  placed 
in  kM  giving  an  excitation  energy  =  (e^.,  +  A^.,)2 

for  the  superconductor.  This  process  can  go  on  only  if 

— * 

the  pair  state  (k'f,  -k’t)  is  initially  empty;  this  occurs 

2 

with  the  probability  u^, .  Energy  is  conserved  if 

|  ^k  |  +  Ek,  =  V.  Another  energy  conserving  process  is 

shown  in  figure  1(b)  ^  where  tunneling  into  k" beneath 

2 

the  Fermi  surface  now  occurs  with  probability  u^n . 

However^  e^-t  =  so  that  from  equations  II-7  and 

2  2 

II-Q,  u^.m  =  Y]^\J  so  ^hat  the  total  probability  for  the 

2  2 

tunneling  process  is  u^,  +  v^,  =  1  as  far  as  Pauli 
principle  restrictions  are  concerned.  The  tunneling 
current  is  then  given  by  summing  only  over  states  above 
(or  below)  the  Fermi  surface  in  the  superconductor  and 
replacing  the  coherence  factor  by  unity. 

The  BCS  theory  has  been  reviewed  in  detail  by  several 
authors,  See3  for  example.,  books  by  J.  R.  Schrieffer 
(1964)  and  G.  Rickayzen  (1965)0 


< 


, 


Figure  1 


Two  final  states  for  a  given  initially 
occupied  state  in  the  tunneling  process 
between  a  normal  and  a  superconducting 
metal.  (Taken  from  Schrieffer  —  "Theory 
of  Superconductivity"  (1964)  p  83.) 
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II .5  Introduction  to  Strong-Coupling  Superconductivity 

The  BCS  theory  contained  no  explicit  assumptions 
about  the  electron-phonon  interaction.  In  fact*  any 
interaction  which  could  produce  the  effect  of  V  would  do 
equally  well  for  their  theory.  This  is  the  only  adjust¬ 
able  parameter  in  the  theory  and  is  to  be  obtained  from 

measured  quantities  such  as  T  ,  etc.  The  laws  of  corres- 

0 

ponding  states  which  arise  naturally  from  this  type  of 
single  parameter  theory  are  satisfied  surprisingly  well 
by  most  superconductors,  so  that  the  form  of  the  inter¬ 
action  is  not  critical  for  many  superconductors.  However, 
for  Pb  and  Hg,  two  superconductors  with  low  values 
and  high  values  of  T  ,  significant  deviations  from  the 
simplest  BCS  predictions  occur. 

The  fact  that  these  deviations  occur  is  not  surpris¬ 
ing  if  one  considers  that  the  weak  coupling  assumption 
kTc  «  ftcDq  is  no  longer  satisfied.  Stating  the  condition 
as  a  ratio  Qjy/T^  we  find  that  for  Al,  a  typical  BCS  super¬ 
conductor,  the  ratio  is  ~350  while  for  Pb,  the  ratio  is 
~l4.  Both  Pb  and  Hg  are  referred  to  as  strong-coupling 
superconductors  and  the  ratio  0^/ Tc  is  often  used  to 
measure  the  amount  of  strong-coupling  character  a  given 
superconductor  has. 

One  can  illustrate  the  deviations  from  BCS  behaviour 
for  strong-coupling  Pb  by  listing  the  following  experi- 


' 

. 


23 


mental  results. 

(1)  The  condensation  energy  at  T  =  0°K  is  given  by 
H^/8tt.  For  Pb  this  equals  24,900  +  erg/cm^0  The 

BCS  expression  for  the  condensation  energy  is  given  as 
yN(0) =  31*000  +  2,000  erg/cm^  (Swihart  et  al  1965) . 

Thus  the  BCS  prediction  is  too  high. 

(2)  The  BCS  energy-gap  transition  temperature  relation 
was  2A  (0)  =  3-528kT  .  In  Pb  the  experimentally  obtained 
value  is  2Aq(0)  =  4.45kTc  and  in  Hg  it  is  2Aq(0)  =  4.6kTc. 
The  value  of  the  gap-T  ratio  shows  in  general  a  strong 
correlation  with  the  C^^/T  ratio,  and  consequently  it 

can  also  be  used  to  measure  the  degree  of  strong-coupling 
in  a  superconductor.  Weak-coupling  superconductors 
generally  show  good  agreement  with  the  BCS  value  of 
2Aq(0)  =  3.53  kTc. 

(3)  The  BCS  theory  predicts  an  electronic  specific  heat 
jump 

ACel  =  1.43  7TC. 

This  gives  ^30.8  mJ  -  mole-'*'  -  °K~^  using  the  measured 
value  of  the  transition  temperature  of  Pb .  Directly 
measured  calorimetric  results  for  the  specific  heat 
jump  vary  from  ^52.9  -  58.1  mJ  -  mole  ^  -  °K  ^  (Clement 
et  al  1952,  Decker  et  al  1958,  and  Shiftman  et  al  1963) . 

(4)  The  critical  field  H  shows  a  positive  deviation 
from  the  simple  relation  =  H  (l  -  t^)  where  t  =  T/Tc 


- 

. 
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is  the  reduced  temperature.  The  3CS  theory  predicts  a 
small  negative  deviation,  which  is  found  in  all  super¬ 
conductors  so  far  investigated  with  the  exception  of  Pb 
and  Hg. 

(5)  The  density  of  states  as  obtained  from  electron¬ 
tunneling  experiments  shows  large  resonance  type  deviations 
from  the  BCS  expression  (eqn.  11-19) . 

II. 6  Zero-Temperature  Strong- Coupling  Theory 

The  BCS  theory  itself  can  be  worked  out  in  the  strong¬ 
coupling  limit  (i.e.  N(0)V  ^  l) .  See,  for  example,  Thouless 
(i960).  It  is  found  that  no  drastic  deviations  from  the 
weak-coupling  limit  occur.  The  gap  ratio  2A  (0)/kTo  goes 
to  a  limit  of  4.00  in  the  unphysical  limit  N(0)V  =  00. 

In  Pb,  N(0)V  ^  0.4  and  for  this  the  ratio  still  remains- 
close  to  3.53. 

Calculations  were  also  carried  out  (Swihart  1962) 
assuming  different  forms  for  the  interaction  V,  represent¬ 
ing  more  closely  the  actual  electron-phonon  interaction. 
These  calculations  were  only  partially  successful.  It  was, 
for  example,  not  possible  to  reproduce  in  detail  the 
observed  structure  in  the  electron  tunneling  data  of  Pb . 

All  these  calculations  use  an  instantaneous  electron- 
electron  interaction,  and  also  neglect  lifetime  effects. 

The  phonon  mediated  electron-electron  interaction  is. 
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however.,  a  retarded  interaction.  In  the  case  of  strong¬ 
coupling.,  this  has  to  be  taken  into  account.  This  leads 
to  qualitatively  different  results  than  those  of  an 
instantaneous  interaction  theory.  One  must  treat  the 
system  using  a  Hamiltonian  which  contains  explicitly  the 
electron-phonon  interaction  (the  so-called  "Frohlich 
Hamiltonian" (Frohlich  1950))  rather  than  representing 
the  interaction  by  an  effective  instantaneous  electron- 
electron  interaction. 

Eliashberg  (i960)  showed  how  the  ground  state  spectrum 
could  be  obtained  for  this  Hamiltonian.  These  calculations 
were  extended,  with  particular  emphasis  on  the  electron¬ 
tunneling  for  Pb,  by  Scalapino,  Schrieffer,  and  Wilkins 
(1963*  1966),  by  Schrieffer  (1964),  and  by  Scalapino  and 
Anderson  (1964).  All  of  these  approaches  use  the  field 
theoretic  Green’s  functions  technique.  The  results  can 
be  expressed  in  terms  of  two  complex  functions  of  energy 
Z(cjd)  and  cp(co)  (where  the  energy  cd  is  measured  relative  to 
the  Fermi  energy).  A  further  function  A(oo)  =  cp(a>)/Z(oo) 
is  introduced.  In  the  weak- coupling  limit  where  this 
theory  goes  over  into  the  weak-coupling  BCS  theory,  A(gd) 
goes  over  into  the  energy  gap  Aq.  In  the  normal  state 
A(gd)  vanishes.  A(cd)  is  therefore  the  analogue  of  the 
energy  gap  in  the  BCS  theory.  There  is  still,  as  in  the 
weak-coupling  theory,  a  certain  minimum  energy  for  creation 
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of  a  quasi-particle *  this  gap  is  denoted  by  A  .  For 
energies  below  the  gap  edge  energy.,  A(cd)  is  constant 
and  therefore*  in  the  strong- coupling  theory*  the  gap 
is  usually  referred  to  as  the  "gap  at  the  gap-edge"*  and 
is  denoted  by  A(A)  =  A  .  Since  A(cd)  is  generally  complex* 
it  is  often  expressed  as  A(cq)  =  A-^(cd)  +  iA^  (  cd)  . 

The  functions  cp(cb)  and  Z(cd)  are  given  by  the 
following  expressions: 


cp(a))  =  /  dcD«  Re ( - 5 - q— i )  (K,  (o)*a)' ) 

o  (cd*2  -  A 1 2 )  2 


Cd[1 


Z(cd)]  = 


c 

/ 


dCD1 


CD  * 

- 0-1  )  K  (CD 

- 


-  N(0)Uc) 
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The  kernels  K+  and  K_  are  given  by: 


00 


K_j_  =  Z  J  dva2(  v)FA(  v)  [  (cd>  +  cd  +  v  +  10+)_1 


A  o 


+  ( CD*  -  cd  +  v  -  iO+)“1].  H-37 


Here  F^(v)  is  the  normalized  phonon  spectrum*  the  polar- 
A 

o 

ization  being  denoted  by  A*  a^(v)  is  an  average  electron- 
phonon  coupling  constant*  cd  is  a  cut-off  energy  taken  at 
approximately  cdc  =  10  cd^*  and  Uc  is  the  Coulomb  pseudo- 
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potential  characterizing  the  short-range.,  instantaneous 
Coulomb  interaction.  The  explicit  appearance  of  the  phonon 
spectrum  and  electron-phonon  coupling  constant  makes  the 
results  of  the  theory  strongly  dependent  on  the  specific 
properties  of  the  system.  Scalapino  et  al  (1965)  have 
calculated  the  s-wave  average  of  the  electron-phonon 

p 

coupling  qC(cd)  .  They  obtain., 


a^(a>) 


=  qD 


.2 
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2  /  dflq 
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q  +  K 
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67 T  K  4tt 


2a> 


qA 


2pw  |  q  +  K 


Q(  2p 
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q  +  k  |)  j  ^ 
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Here  q^  is  the  Debye  wave  vector.,  pp  the  Fermi  wave  vector, 
and  m  is  the  mass  due  to  the  band  structure  with  Coulomb 
corrections  but  without  corrections  for  the  electron- 
phonon  interaction.  9  is  a  unit  step  function  vanishing 
when  its  argument  is  less  than  zero  and  equal  to  1  other¬ 
wise.  A  discussion  of  the  results  of  the  theory .,  must 
therefore  be  based  on  specific  assumptions  regarding  the 
system  under  consideration. 

The  strong-coupling  theory  has  also  been  developed 
for  finite  temperatures  using  thermodynamical  Greens 
functions  techniques.  Some  of  the  results  will  be  given 
in  the  next  sections. 


. 
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I I . 7  Application  to  the  Tunneling  Process 

Scalapino„  Schrieffer„  and  Wilkins  (1966)  and 
Schrieffer,  Scalapino„  and  Wilkins  (1963)  used  strong¬ 
coupling  theory  to  calculate  the  tunnel  current  in  a 
superconductor-normal  metal  tunnel  junction.  They  found 
that  equation  11-33  still  holds  but  now  the  gap-parameter 
is  replaced  by  the  energy  dependent  gap  A (00)  such  that 


03 

(CD2  -  A(C0)2A 

Thus  strong-coupling  theory  replaces  the  constant  BCS 
energy  gap  by  the  function  A (03)  and  the  normalized 
tunneling  density  of  states  becomes  Nrp(a3)/N(0)  .  This  is 
in  1:1  correspondence  with  the  experimentally  measured 
junction  conductance  as  before.  The  normalized  density 
of  states  can  be  calculated  by  solving  simultaneously 
the  coupled  equations  for  A(cd)  for  a  particular  choice 
of  phonon  spectrum  and  coupling  constant  a  (cjd)  and  then 
using  the  energy  dependent  value  of  A (03)  in  equation  11-39* 
Calculations  along  these  lines  have  been  performed  by 
Rowell  and  Kopf  (1965);  by  Scalapino  and  Anderson  (1964),, 
and  by  others.  In  general  these  calculations  showed 
excellent  agreement  with  experimentally  obtained  curves. 

McMillan  and  Rowell  (1965)  took  the  further  step  of 

2 

actually  inverting  this  process.  They  guess  at  a  (cd)f(o)) 


II-39 


NT(a3) 

— -  =  Re  ( 

N(0) 
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and  solve  the  gap  equations  choosing  U  to  give  the 
experimental  value  of  A  .  They  then  calculate  the 
electronic  density  of  states.  Linear  feedback  techniques 
are  then  used  to  correct  a  (cd)F(cd)  for  the  "error" 

Naxp(cu)  -  N^a"1'C(a))  .  In  this  way  they  can  obtain  the 

O  O 

actual  form  of  a  (a>)F(a>)  from  experimental  data  on  the 
density  of  states.  Their  result  Is  shown  In  figure  2 
along  with  plots  of  the  normalized  conductance  a  (i.e. 
the  density  of  states)  and  da/dV. 

II. 8  Phonon  Emission  Processes  and  the  Density  of  States 

The  structure  in  the  normalized  conductance  curve 
(as  well  as  in  da/dV)  centered  around  4.4  meV  and  8.5  meV 
is  associated  with  phonon  emission  processes.  Physically 
these  processes  occur  in  the  following  way.  Quasi-particles 
at  the  longitudinal  and  transverse  energies  co^  and 
have  a  relatively  high  chance  of  losing  energy  by  emission 
of  phonons.  They  then  drop  to  the  peaked  density  of 
states  region  near  the  gap  edge.  Energetically  this 
occurs  when 

V  =  A  +  cd,  11-40 

O  A 

and  to  a  lesser  extent  at  harmonic  frequencies. 

The  resonance  structure  is  reproduced  in  model  cal¬ 
culations  for  the  density  of  states.  This  can  be  seen 


- 
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Figure  2 


The  normalized  conductance  Oj  and  da/dV 
for  a  Pb-Al  tunnel  junction.  Also  shown 
is  a  (<d)F(cd)  calculated  by  McMillan  and 


Rowell  (1965). 
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most  easily  by  expressing  A(co)  as 

A  (go)  =  A1(cd)  +  iA2(co)  II-41 

where  A^  (co)  and  Aq(od)  are  real  functions.  These  functions 

are  then  calculated  for  a  model  chosen  to  represent  Pb. 

Figure  3  shows  curves  of  A^fco)  and  A2(co)  for  T  =  0  and  a 

particular  model.  If  we  now  expand  the  density  of  states 

2 

to  first  order  in  A  we  have 

Nfo))  A^(cd)  -  Ap(oj) 

— -  =  1  +  — - — — - +  ,  .  .  11-42 

N(0)  2cd 

The  emission  process  is  reflected  as  structure  in  the 
plots  of  A^fco)  and  Ap(co),  and  from  ean.  11-42,  in  the 
density  of  states.  Just  below  the  transverse  phonon 
emission  threshold  A2(cd)  increases  writh  increasing  co, 
causing  Nrp(co)  to  rise.  After  the  threshold,  A^  decreases 
while  A2(co)  becomes  large,  thereby  producing  a  sharp  drop- 
in  Nrp(co)  .  A  similar  situation  occurs  near  the  longitu¬ 
dinal  phonon  emission  threshold.  Thus  a  good  deal  of 
detailed  information  about  the  phonon  density  of  states 
can  be  obtained  from  tunneling  curves.  In  particular, 
a  study  of  figure  2  shows  that  the  peaks  in  a  (cd)f(go) 

(that  is,  peaks  in  f(cd))  occur  at  the  same  energies  as 
minima  in  dcr/dV,  and  occur  at  the  mid-height  points  of 
the  structure  drop-offs  in  a.  One  should  therefore  be 
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Figure  3 


A  plot  of  A-^(cd)  and  A^cd)  for  reduced, 
temperatures  of  t  =  0  and  t  =  0.9&. 
(Taken  from  Scalapino,  Wada^  and  Swihart 

(1965) .) 
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able  to  use  the  conductance  structure  and  da/dV  for 
measuring  phonon  peak  shifts  under  pressure.  This  is 
discussed  further  when  such  results  are  presented  In 
section  VI. 2. 

McMillan  and  Rowell  (1965)  have  used  their  calculated 
results  of  a  (a))F(a>)  In  estimating  the  effective  mass 
of  Pb.  This  was  done  by  calculating  the  real  part  of 
the  renormalization  parameter  Z(cd)  for  a)  =  0  and  T  =  0. 
They  then  have  (since  A  =  0  as  well), 

p 

m  co  or  (as) F( cd) 

Z  (0)  =  —  =1  +  2  / - den 

m  oo) 

■X 

where  m  is  the  apparent  mass  due  to  electron-phonon 
renormalization  o  Using  their  plot  of  ol~  (cd)F(o))  they 
obtain 

m  =  2 .33  +  0 .02  m 

where  m  is  the  ordinary  electron  mass.  This  compares 

with  2.2  and  2.4  obtained  theoretically  by  Swihart, 

Scalapino,  and  Wada  (1965)^  using  model  calculations  of 
o 

a  (a))F(a>)  for  Pb,  and  with  a  value  of  2.1  obtained  by 
specific  heat  measurements  (Keesom  et  al,  1963) .  The 
value  of  the  effective  mass  (which  can  in  large  measure 
be  accounted  for  by  the  electron-phonon  interaction) 
demonstrates,  in  yet  one  more  way,  the  importance  of 
the  strong  interaction  on  the  properties  of  Pb. 
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II. 9  Strong-Coupling  Theory  at  Finite  Temperature 

Scalapino,  Wada,  and  Swihart  (1965)  and  Swihart, 
Scalapino,  and  Wada  (1965)  have  extended  the  gap  calcu¬ 
lations  to  finite  temperature  using  strong-coupling 
theory.  They  have  published  curves  of  A^(u>,  T)  and 
A^cq,  T)  for  t  =  T/Tc  =  O.98  and  t  =  0.95  and  a  model 
chosen  to  represent  Pb.  In  these  calculations  they 
chose  for  Pb  a  lattice  spectrum  consisting  of  two 

Lorenztians  centered  at  o^_  =  4.4  meV  and  cjd  =  8.5  meV. 

2  2 

The  electron-phonon  coupling  constants  and  were 

chosen  to  obtain  agreement  with  the  experimentally 

determined  zero  temperature  gap  at  the  gap-edge. 

2  2 

Calculations  were  made  for  a^/a^  =  1  (case  1)  and 
arat  =  ^  (case  2).  The  t  =  0 . 98  results  are  included 
in  figure  3*  Their  results  can  be  summarized  as  follows 

(1)  The  resonance  peaks  connected  with  phonon  emission 
are  sharpened  at  non-zero  temperatures,  and  are  shifted 
down  in  energy  in  accordance  with  a  reduction  in  the 
energy  gap  with  increase  in  temperature . 

(2)  Near  cd  =  0,  A-,  (cjo)  is  flatter  as  a  function  of 
energy. 

(3)  New  resonance  structure  occurs  near  the  frequencies 
ox^  -  Aq(t)  which  is  most  evident  near  u,_  -  Aq(t), 
where  A2(o), T)  goes  through  a  negative  resonance  peak. 


. 
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The  new  resonances  are  thought  to  be  associated  with 

recombination  processes  between  injected  electrons  (as 

from  tunneling)  and  thermally  excited  quasi-particles 

in  the  high  density  of  states  region  near  the  gap  edge 

A  (T)  .  The  combined  particles  emit  a  phonon,,  and  form 

a  ground  state  pair.  The  rate  of  the  process  is  a 

maximum  when  the  total  energy  V  of  the  injected  electron 

and  excited  quasi-particle  at  A  (T)  is  equal  to  a 

frequency  corresponding  to  a  peak  in  the  phonon  density 

of  states  (i.e.  V  +  A  (T)  =  ) .  This  additional  struc- 

o  A 

-  ture  in  A(oj)  should  appear  in  the  tunneling  density  of 
states.  They  have  published  theoretically  expected 
tunnel  junction  conductance  curves  for  t  =  0.95  and 
t  =  0.92.  An  investigation  of  this  effect  is  discussed 
in  section  VI. 6. 

(4)  The  calculation  of  Tq  from  the  gap  relation  leads 
to  2A  (0)  =  4 .40  kT,  in  case  1  and  2A  (0)  =  4.33  kT 
in  case  2. 

(5)  The  reduced  gap  at  the  gap-edge  energy  A  (T)/A  (0) 
as  a  function  of  reduced  temperature  t  is  "within  a 
couple  of  percent"  of  the  BCS  expression. 

An  additional  comment  on  the  gap-T  relation  is  in 
order.  Using  the  measured  gap  value  for  Pb  (A  (0)  =  1.34 
meV) }  the  transition  temperature  would  have  to  be  ^  11.3°K 
for  the  BCS  ratio  2A  (0)/kT  =  3.53  to  hold.  Wada  (1964a) 
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associates  the  observed  value  of  ^  4.4  with  the  strong¬ 
coupling  resonant  processes  which  give  rise  to  quasi¬ 
particle  damping.  Since  the  damping  rate  is  greater 
at  higher  temperatures  where  more  excited  quasi-particles 
are  produced*  the  transition  temperature  is  reduced  much 
more  than  the  zero  temperature  energy  gap*  thereby 
increasing  the  ratio  2AQ(0)/kTc. 

II. 10  Calculation  of  Thermodynamic  Quantities 

Wada  (  1964b)  has  undertaken  a  .calculation  of  some 
of  the  thermodynamic  quantities  in  Pb  using  strong¬ 
coupling  theory.  In  particular*  he  attempts  to  reconcile 
the  discrepancies  in  AC  ^  and  the  condensation  energy 
relation  mentioned  earlier.  He  has  shown  that  in  general 


N(0)  I(T) 


11-43 


8tt  4a  -  1 

where  H  is  the  zero-temperature  critical  field*  and 


o 


l(T)  is  a  functional  of  the  electron-phonon  renormal¬ 
ization  functions  and  Zg  and  of  the  complex  gap. 
The  exponent  a  appears  in  the  isotope  effect  as 


H  oc  M  a 
o 

and  a  =  \  for  Pb 0  Therefore*  4a  -  l  reduces  to  1 
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so  that 


H 

—  =  N(0)I(T) 
87T 
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In  the  weak  coupling  limit*  I  goes  properly  over 

i  2 

into  I  =  y Aq  to  yield  the  BCS  result.  One  can  there- 

fore  introduce  the  ratio  I/J-Ag  and  use  this  as  an 

approximate  measure  of  the  coupling  strength.  For  Pb* 

2 

l/§-A^  =  0.83  so  the  condensation  energy  is  smaller  than 
the  BCS  prediction. 

Wada  has  also  derived  an  expression  for  the  jump 
in  electronic  specific  heat: 


AC  =  -  N(0)  d  I(T) 
el  '  - p — 

t-1  dr 
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Carrying  out  the  derivative  of  the  function  I  he  obtains* 
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where  p  =  1/kT,  A  =  A(<d),  and  k  is  Boltzmann's  constant. 
Thus  knowing  the  proper  energy  dependence  of  8a  /86 

0 

one  can  in  principle  calculate  AC  .  Wada  remarks*  "we 
may  conclude  that  the  present  theory  is  likely  to  predict 
a  larger  jump  in  specific  heat  and  give  better  agreement 
with  the  experimental  results  for  the  strong  coupling 
superconductors  than  the  BCS  results." 
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CHAPTER  III 

THE  CRYOSTAT  AND  PRESSURE  SYSTEM 

III.l  The  Cryostat 

The  cryostat  was  suspended  inside  a  two-dewar 

coolant  system.  The  outer  dewar  contained  liquid 

nitrogen  and  the  inner  dewar  contained  liquid  helium 
z  4\  i 

four  (He  ) .  Figure  4  shows  a  cross-section  of  the 
cryostat . 

Since  a  working  temperature  range  of  '■'ol  •  4°K  to 
i  o  4 

a,40  K  was  required,  a  He  reservoir  was  soldered  to  the 
specimen  bomb.  By  admitting  liquid  helium  from  the 
inner  dewar  through  the  needle  valve  control  and  then 
reducing  pressure  on  the  enclosed  liquid  vapor  with  a 
vacuum  pump,  a  low  temperature  limit  of  approximately 
1.25°K  could  be  reached.  Installation  of  heater  resis¬ 
tors  enabled  the  higher  temperatures  to  be  reached 
quickly  when  necessary. 

The  cryostat  vacuum  pump  tubes  consisted  of  3/^ 
inch  diameter  thin-wall  stainless  steel  tubing.  This 
was  chosen  because  of  its  strength  and  low  thermal  con¬ 
ductivity.  The  high  pressure  capillary  tube  entered 
the  cryostat  vacuum  jacket  ten  inches  above  the  top  of 


* 


. 


Figure  4 


A  cross-section  of  the  cryostat  used0 
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the  vacuum  can.  '  By  inserting  stainless  steel  tubing 
between  the  vacuum  can  and  the  capillary  contact,  the 
pressure  tubing  was  thermally  isolated  from  the  vacuum 
can.  This  step  was  taken  to  help  prevent  capillary 
blockage  at  the  point  of  contact  during  pressurization. 

Manganin  resistance  wire  was  wound  spirally  around 
the  outside  of  the  pressure  capillary  both  inside  and 
outside  the  vacuum  can.  By  passing  a  current  through 
these  wires,  extra  heat  was  supplied  to  the  capillary 
to  help  prevent  its  blockage 0  Even  after  these  measures 
were  taken,  freezing  of  the  high  pressure  gas  at  the 
capillary  inlet  to  the  vacuum  can  still  posed  a  definite 
problem . 

4  / 

The  He  melting  curve  (figure  5)  shows  that  at  a 
pressure  of  5^000  atmospheres,  the  freezing  temperature 
T  ,  is  quite  high  at  ^40°K.  After  pressurizing  the  bomb, 
about  600  cc  of  liquid  helium  from  the  inner  dewar  was 
needed  for  freezing  this  gas  at  constant  pressure. 

During  the  earlier  runs  it  was  found  that  if  this  liquid 
was  in  contact  with  the  vacuum  can,  the  high  pressure 
capillary  invariably  blocked. 

The  problem  was  solved  by  connecting  stainless  steel 
tubing  to  the  needle  valve  control  leading  to  the  helium 
reservoir.  By  extending  this  tubing  about  eight  to  ten 
inches  below  the  base  of  the  vacuum  can,  which  was  now 
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Figure  5 


4  / 

The  He  melting  curve.  (Due  to  Mills 

and  Grilly  (1955)  and  Dugdale  (1965)0) 
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no  longer  In  contact  with  the  liquid  helium.,  sufficient 
helium  coolant  could  be  sucked  into  the  reservoir  to 
produce  the  freezing.  The  capillary  now  remained  open 
so  constant  pressure  freezing  occurred.  This  resulted 
in  frozen-in  pressures  about  twice  as  large  as  constant 
volume  freezing  could  produce.  The  pressure  drop  between 
T^  and  the  working  temperature  was  about  10$  which  is  in 
good  agreement  with  published  isochores  (figure  6) . 

Thermometer  leads  were  fed  into  the  vacuum  can 
through  a  glass  "cow"  vacuum  seal  provided  with  kovar 
feed-throughs.  This  method  insured  adequate  lead  anchor¬ 
ing  at  4.2°K  since  the  leads  passed  through  several 
inches  of  liquid  helium  in  the  inner  dewar .  Leads  for 
the  capillary  heater  were  fed  down  one  of  the  vacuum 
pump  tubes  and  anchored  to  a  copper  post  soldered  to 
the  vacuum  can. 

III. 2  The  High  Pressure  Bomb 

The  bomb  used  for  the  experimental  pressure  chamber 
is  illustrated  in  figure  7.  The  bomb  material  is  copper 
with  about  2$  beryllium  added  to  increase  its  hardness. 
This  material  is  preferrable  to  stainless  steel  because,, 
while  both  have  the  same  approximate  rupture  strength, 
the  Cu-Be  has  a  much  higher  thermal  conductivity,  and  does 
not  become  brittle  at  low  temperature.  The  Cu-Be  was 
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Figure  6 


He  isochores.  (Taken  from  Dugdale 


and  Simon  (1953).) 
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Figure  7 


The  high  pressure  bomb. 
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purchased  as  a  half  hardened  alloy  from  "The  Beryllium 

Corporation",  Reading,  Pennsylvania.  After  machining.  It 

was  fully  hardened  by  heating  it  to  a  temperature  of  600°K 

for  l\  hours.  A  diamond  press  hardness  test  of  the  full 

and  half  hardened  samples  verified  that  the  hardness  was 

indeed  about  twice  as  great  after  the  heat  treatment. 

The  lense  ring  which  forms  the  pressure  seal  between 

the  nut  and  the  bomb  was  made  of  stainless  steel,  ground 

and  polished  with  edges  cut  at  37°.  Hence  an  increase 

in  gas  pressure  increases  the  outward  force  on  the  ring 

so  that  it  acts  as  a  self  tightening  seal.  Notice  the 

relief  hole  drilled  into  the  bomb  to  release  gas  that 

might  escape  past  the  lense  ring. 

4 

The  cylindrical  He  reservoir  was  soldered  around 

the  bomb.  The  inlet  tube  from  the  needle  valve  goes  to 

the  bottom  of  the  reservoir.  Thus  during  the  first  stages 
4 

of  He  admission,  the  bottom  of  the  bomb  was  the  coldest. 
Therefore,  the  solid  helium  formed  from  the  bottom  of 
the  bomb,  and  then  grew  toward  the  top. 

Thermal  isolation  between  the  pumping  chamber  and 
the  vacuum  can  was  achieved  by  inserting  about  10  inches 
of  spiral  cupro-nickel  into  the  lines  joining  the  needle 
valve  and  the  vacuum  pump  tube  to  the  reservoir. 

A  heater  resistor  was  attached  to  the  top  of  the 


. 
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bomb  nut  and  a  calibrated  germanium  thermometer  was 
mounted  on  the  reservoir.  Thermometer  lead  anchoring  to 
the  specimen  bomb  was  achieved  by  wrapping  the  leads 
around  the  reservoir  several  times  and  then  cementing 
them  to  it  with  General  Electric  baking  varnish  type 
GE9825 .  The  lowest  temperature  attainable  in  the  cryostat 
was  about  1.25°K. 

The  specimen  leads  were  introduced  into  the  bomb 
through  the  high  pressure  gas  capillary.  The  leads  vjere 
terminated  on  a  teflon  strip  as  shown.  Specimen  mounting 
involved  soldering  the  specimen  leads  to  the  appropriate 
terminals . 

III. 3  Pressure  Production 

The  gas  pressure  system  was  purchased  from  Harwood 
Engineering.,  Walpole,  Massachusetts .  It  consists  of  three 
pumping  stages.  The  first  stage  and  the  second  stage 
occur  at  either  ends  of  a  double  acting  piston  which  is 
oil  driven  and  reversed  for  consecutive  strokes  by  elec¬ 
trical  micro-switches  which  control  an  oil  diverting 
mechanism.  Stage  1  exhausts  into  stage  2.  An  inlet 
pressure  of  800  psi  at  stage  1  is  increased  to  a  maximum 
of  12,000  psi  at  the  output  of  stage  2.  This  stage 
exhausts  into  an  intensifier  chamber  which  can  be 
isolated  from  the  stage  2  pump  by  manually  operating  a 
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needle  valve.  The  third  or  intensifier  stage  Is  an  oil 
driven  piston  stage  (area  multiplication  of  14) .  The 
oil  pressure  is  provided  by  an  air  diaphragm  pump  for 
this  stage.  A  motor  driven  pump  provides  oil  pressure 
to  stages  one  and  two.  The  intensifier  is  connected 
through  another  needle  valve  to  the  specimen  bomb. 

The  bomb  pressure  (and  intensifier  pressure  if  the 
connecting  needle  valve  is  open)  is  measured  by  using 
the  resistance  of  a  manganin  wire  as  the  sensing  device. 
The  resistance  is  measured  with  a  Wheatstone  bridge 
circuit  and  the  imbalance  is  calibrated  in  units  of  psi. 
The  instrument  accuracy  is  better  than  1%. 

III. 4  High  Pressure  Electrical  Lead  Seal 

Six  electrical  leads  to  the  inside  of  the  high 
pressure  bomb  were  needed.  This  allows  four  terminal 
measurements  to  be  made  on  two  junctions  mounted  on  the 
same  substrate.  The  high  pressure  lead  seal  described 
by  Goree,  McDowell,  and  Scott  (1965)  was  used.  This 
method  involves  a  room  temperature  high  pressure  seal.  . 


Size  34  gauge  insulated 
copper  wires  were  fed 
through  about  12  feet  of 


high  pressure  capillary  and 


couplings  as  illustrated. 
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Eccobond  104  epoxy  purchased  from  Emerson  and  Cumming 
Inc.  was  squeezed  into  the  high  pressure  tube  at  the 
two  places  shown.  Before  pressurizing,,  oil  is  admitted 
and  frozen  in  the  U-tube  with  liquid  .  We  have  found 
that  in  this  way  pressure  up  to  60,000  psi  could  be 
contained  quite  adequately.  However,  after  two  or 
three  months  of  standing  at  room  temperature  and  zero 
pressure,  the  epoxy  seals  had  a  tendency  to  soften  and 
fail  under  subsequent  pressurizing.  We  ascribe  this  to 
the  presence  of  the  oil.  For  pressures  near  100,000 
psi,  it  would  probably  be  necessary  to  make  a  new  seal 
before  every  run.  For  pressures  around  50,000  psi,  a 
plug  of  frozen  oil  would  hold  the  pressure  even  if  the 
gas  leaked  past  the  epoxy  seals. 

III. 5  Precautionary  Measures 

The  production  of  high  gas  pressure . coupled  with 
low  temperature  requires  some  safety  precautions  since 
a  vacuum  failure  can  result  in  a  rapid  warming  up  of  the 
system  with  a  consequent  large  pressure  build-up.  In 
assembling  the  pressure  system  the  following  precautions 
were  taken; 

(l)  Pressure  lines  were  either  shielded  by  angle 
iron  or  tied  down  to  prevent  possible  whip-lashing  if 
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one  of  the  joints  broke. 

(2)  All  joints  were  positioned  so  that  they 
pointed  away  from  potential  working  areas . 

(3)  All  capillary  connection  joints  were  conical 


A  vent  between  the  seal  and  thread  was  always 
present  to  prevent  pressure  build-up  in  the  event  of  a 
leak  at  the  cone  seal.  Absence  of  this  vent  could  have 
resulted  in  ejection  of  the  tube  if  a  leak  at  the  seal 
developed . 

(4)  The  whole  cryostat  was  mounted  inside  a  1/4 
inch  plate  steel  closet  lined  inside  and  out  with  3/4 
inch  plywood.  This  completed  the  shielding  of  all  of 
the  high  pressure  volume. 

III. 6  Thermometry 

Temperature  measurements  were  obtained  using  a 

*  . 
germanium  thermometer  .  This  was  calibrated  to  0.1% 

* 

I  wish  to  thank  Dr.  Franck  for  the  use  of  his  cali¬ 
brated  thermometer  throughout  this  work. 
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3  4 

accuracy  against  He  and  He  vapor  pressures  from 
0.65°K  to  4.2°K.  For  temperatures  from  4.2°K  to 
~7*5°K  the  thermometer  was  calibrated  against  a  second 
germanium  thermometer,  which  in  turn  had  been  cali¬ 
brated  against  a  gas  thermometer  to  an  accuracy  of  0.1 % . 
During  the  calibration  of  the  second  germanium  ther¬ 
mometer,  a  re-determination  of  the  transition  tempera¬ 
ture  of  unstrained  pure  bulk  Pb  gave  J. 196  +  0.005°K, 
in  good  agreement  with  the  value  of  7*193  +  0.005°K 
reported  by  Franck  and  Martin  (1961) .  The  estimated 
temperature  uncertainty  is  +  0.1 %  for  all  tempera¬ 
tures  below  7-5°K. 
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CHAPTER  IV 

THE  ELECTRONIC  DETECTION  SYSTEM 

IV . 1  Introduction 

The  investigation  of  characteristic  microscopic 
metal  properties  at  low  temperature  calls  for  especially 
low  level  probe  signals.  Ideally,  thermal  smearing,  kT, 
should  be  the  limitation  where  k  is  Boltzmann's  constant. 
Since  kT  is  approximately  86  peV  at  1°K,  our  boundary 
condition  is  clearly  defined. 

Noise  arises  from  such  statistical  fluctuation 
phenomena  as  Johnson  noise  in  resistors  and  shot  noise 
in  vacuum  tubes  and  semiconductors;  these  both  produce 
a  white  noise  spectrum  in  which  the  noise  power  per 
unit  bandwidth  is  the  same  at  all  frequencies .  Another 
source  is  the  gain  modulation  or  flicker  effect  noise 
associated  with  both  vacuum  tubes  and  amplifying  circuits. 
This  noise  frequency  spectrum  varies  as  l/f  so  that  a 
large  contribution  to  the  total  noise  occurs  near  dc . 
Interference  phenomena  which  are  not  really  noise,  but 
produce  the  same  effect  in  obscuring  the  signal,  include 
power  line  pick-up  and  radio-frequency  interference. 

With  respect  to  white  noise,  little  can  be  gained 
by  moving  the  signal  frequency  to  a  different  value. 
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However.,  the  contribution  of  the  white  noise  to  the 
output  voltage  of  the  detection  system  is  directly 
proportional  to  the  square  root  of  the  bandwidth,  and 
can  be  reduced  to  an  arbitrarily  small  value  by  reducing 
the  bandwidth.  For  l/f  noise  and  interference,  an 
operating  frequency  different  from  dc  and  from  inter¬ 
fering  frequencies  can  be  selected.  Considerable  effort 
was  made  to  shield  against  pick-up.  The  residual  level 
was  about  1  to  2  microvolts, . and  consisted  of  60  cycle 
and  1  megacycle  contributions.  The  latter  contribution 
was  the  smaller  and  came  from  the  computer  stationed 
four  floors  away. 

The  tunnel  junction  normalized  conductance  a  is 
related  to  the  superconducting  density  of  states  by  the 
relation  a  =  (di/dv)ns/(di/dv)nri  =  Nt(<d)/N(0)  .  Both  a 
and  its  derivative  da/dv  can  be  obtained  by  applying  an 
ac  signal  and  using  harmonic  detection  techniques. 

If  the  voltage  modulation  amplitude  is  kept  constant, 
the  current  i,  passing  through  the  junction  may  be  written 
in  terms  of  a  Taylor  series: 

i(v)  =  i(v  )  +  (di/dv)  m  cosaut  +  -^(d^i/dv^)v  m2(l  +  cos2ayt) 

o  o 

.  IV- 1 

T . . 

where  m  is  the  constant  voltage  modulation  amplitude,  and 
vq  is  the  dc  junction  bias  voltage.  Therefore,  the  com-* 
ponent  of  voltage  across  the  junction  at  angular  fre- 
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quency  co  is  proportional  to  di/dv  and  the  component  at 

2  2 

2o)  is  proportional  to  d  i/dv_. 

The  respective  informations  at  frequencies  <d  and  2oo 
were  obtained  using  a  bridge  system  (to  be  described) 
and  a  lock-in  amplifier  as  the  signal  detector. 

A  lock-in  amplifier  is  essentially  a  narrow  band 
detector  in  which  a  signal  is  beat  with  a  reference  signal 
of  the  same  frequency  giving  a  dc  output;  the  upper  side¬ 
band  is  stopped  by  a  low  pass  filter.  The  low  pass 
filter  determines  the  bandwidth  of  the  lock-in  output. 

The  main  advantages  of  the  lock-in  device  over  a  tuned 
amplifier  circuit  are: 

(1)  it  is  easily  tunable  since  the  reference  frequency 
determines  the  center  frequency  of  the  pass  band 

(2)  no  matter  how  narrow  the  bandwidth  of  the  detec¬ 
tion  system,,  the  center  of  the  pass  band  is  always 
locked  to  the  signal  frequency  when  the  signal  is  avail¬ 
able  for  use  as  the  reference. 

Thus  by  moving  the  signal  frequency  to  about  1700 
cycles  per  second.,  and  using  the  lock-in  amplifier  as  a 
detector  with  an  input  frequency  pass  band  of  about  100 
cps,  the  thermal  noise  voltage  is  almost  completely  re¬ 
moved  from  the  signal  input  voltage.  Only  the  thermal 
noise  associated  with  the  pass  band  frequencies  contri¬ 
butes  to  the  noise  background . 
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IV. 2  Circuit  Operation 

A  block  diagram  of  the  detection  system  is  shown  in 
figure  8.  For  the  conductance  measurements,,  signal  at 
frequency  gd  is  required.  In  the  case  of  da/dv,  the  sec¬ 
ond  harmonic  component,  2a>,  generated  across  the  non¬ 
linear  tunnel  resistance  R,.,  must  be  measured.  The  latter 
is  accomplished  by  setting  the  reference  frequency  to  2a> 
in  a  frequency  doubler  circuit  connected  to  the  oscillator. 
At  the  same  time,  the  lock-in  input  frequency  is  set  to  2co. 

The  bridge  circuit  which  produces  the  required  con- 

* 

start  voltage  condition  is  shown  in  figure  9.  The  junction 
resistance  and  capacitance  are  balanced  by  and  R^. 

The  latter  consists  of  two  10  turn  potentiometers  in 
series.  As  the  dc  bias  is  swept  with  a  motor  drive, 
changes  in  the  tunnel  characteristic  results  in  a  voltage 
difference  being  fed  to  the  input  of  the  operational 
amplifier.  This  amplifier  has  a  gain  of  approximately 
50,000.  The  voltage  difference  is  amplified  and  current 
fedback  through  the  10k  feedback  resistor  to  the  junc¬ 
tion.  The  feedback  current  is  such  that  the  signals 
across  R,  and  R  are  equalized  to  a  good  approximation. 


It  is  a  pleasure  to  thank  Dr.J.S.  Rogers  for  his  sug¬ 
gestions  regarding  the  electronics,  and  particularly  for 
his  construction  of  this  bridge  circuit. 
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Figure  8 


A  block  diagram  of  the  detection  system „ 


Figure  9 


The  bridge  circuit. 
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The  feedback  current  required  for  equalization  is  a 
direct  measure  of  di/dv*  as  shown  belowQ  This  is  there¬ 
fore  converted  into  a  voltage  and  fed  to  the  input  of 
the  lock-in  amplifier. 

IV . 3  Bridge  Circuit  Analysis 

Consider  for  the  purpose  of  analysis  the  simplified 
circuit  of  figure  10  (the  deleted  resistors  simply  act 
as  Isolation  resistors  and  current  limiters) .  Assume 
that  =  C^_  so  the  reactance  is  zero.  ( ^  0.025  Rf 
and  thus  dominates  stray  inductance.)  The  current  i 
consists  of  dc  as  well  as  ac  contributions ,  and  the  oper¬ 
ational  amplifier  input  is 


eo  =  A(lRd  -  lRt  - 

61Rfc)  . 

IV- 2 

Also 

e0  =  A(e+  -  e_) 

IV- 3 

=  61 (Rf  +  Rfc)  + 

lRt. 

IV- 4 

Since  the  input  impedance  of  the  operational  amplifier 
is  ^  100  kO.,  and  is  therefore  much  greater  than  no 
circuit  loading  occurs.  Therefore., 

A(i(Rd  -  Rt)  -  5iRt)  =  6i(Rf  +  Rfc)  +  iRt .  IV- 5 

This  reduces  to 

i[Rt(-l  -  A)  +  ARd]  . 

Rf  +  Rt  +  ARt 
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Figure  10 


The  simplified  bridge  circuit 0 
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Setting  P  equal  to 

FL 


t 


P  = 


one  has 


Rf  +  Rt 


IV- 7 


6i  = 


i[6(-l  -  A)  +  PARd/Rt] 


1  +  PA 

iPA  -1  R 

-  [  _  _  1  +  _H] 

1  +  PA  A  R, 


d 

t 


iPA 


1  +  PA 


[Rdgt  -  1], 


iv- 8 


iv- 9 


IV- 10 


where  the  tunnel  conductance  g,  =  l/R,  and  A  ^  50.,  000. 

U  1/ 

The  value  of  PA  varies  from  ^250  to  ^10,000  with  the 

circuit  elements  used.,  and  the  variation  of  R  .  Therefore., 

PA 

to  a  good  approximation  d  +  ^  1,  and  one  then  has 


6i  'v  i[Rdgt  -  1]  . 

Going  back  to  expression  IV-4  for  the  output  voltage  e 
and  substituting  the  value  for  51,  gives 


IV- 11 


o 


eo  -  1[RdSt  -  ^f  +  +  gt 

=  iRdRfgt  +  i(RH  -  R^) . 
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IV- 12 
IV- 13 


Therefore  eQ  is  linear  in  the  junction  conductance  g^,  and 
the  normalized  conductance  a  is  obtained  from 

St 


sn 


a  = 


IV- 14 
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Reactive  and  resistive  balance  of  the  bridge 
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circuit  can  be  sensitively  achieved  using  this  detection 
system.  The  lock-in  device  provides  continuous  (0°  -  360°) 
phase  shift  between  the  reference  signal  and  the  input 
signal.  The  circuit  is  completely  balanced  for  only  two 
values  of  the  phase  setting  (angle  and  supplementary 
angle) .  At  these  settings  a  large  reactive  imbalance 
does  not  affect  the  resistive  balance  and  vice  versa. 

Any  other  phase  setting  results  in  a  false  balance  so  that 
a  change  in  either  the  reactive  or  resistive  value  now 
results  in  a  shift  in  both  balance  readings. 

IV. 4  Conductance  Calibration 

The  ac  signal  voltage  across  R^_  was  usually  set  at 
^  150  pvolts  peak  to  peak  (pp)  or  less.  This  corresponds 
to  -v  50  pvolts  rmSj  or  less  than  1°K  thermal  smearing. 
Modulation  levels  as.  low  as  50  p/volts  pp  were  often 
used . 

The  conductance  was  calibrated  at  a  junction  bias 
of  30  millivolts  because  it  was  found  that  in  this  way 
the  calibration  remained  the  same  whether  the  Pb  was 
normal  or  superconducting.  At  lower  biases^  possible 
variations  in  the  density  of  states  for  the  supercon¬ 
ducting  Pb  might  affect  the  calibration  and  thus  invali¬ 


date  the  results. 
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In  calibrating  the  conductance,  R^  and  were 
first  set  to  give  a  complete  balance 0  The  junction  was 
then  open  circuited  with  the  switch  (figure  9)  in  the 
"zero"  position.  The  signal  from  R^  could  no  longer  be 
equalized  and  the  lock-in  output  now  corresponded  to 
zero  conductance.  The  difference  between  g^  =  0  and 
g^_  =  1  was  then  adjusted  to  give  a  specific  Y  deflection 
on  the  X-Y  recorder.  By  changing  the  lock-in  and  Y-axis 
sensitivities,  fractional  percentages  of  unity  conduct¬ 
ance  could  be  dialed  directly. 

Linearity  tests  on  the  conductance  were  carried  out 
with  the  four  terminal  network  shown  below. 


Lead  resistances  of  10  ohms  were  chosen  to  approximate 
actual  values  o  R^  was  varied  from  100  -  5 .>0000  with 
various  accompanying  values  of  .  The  bridge  was  found 
to  be  linear  in  conductance  to  within  the  Ink  width  of 
the  plotter  pen.  We  estimate  that  the  conductance 
inaccuracy  is  less  than  0.2^  for  all  of  these  values „ 

For  values  of  R^  below  about  1000  the  junction  capaci¬ 
tance  began  to  deviate  from  in  the  reactive  bal¬ 
ance.  This  was  found  to  be  due  to  the  presence  of  the 
lead  resistors.  However,  since  the  linearity  of  the 
conductors  was  of  primary  importance,  this  did  not 
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pose  any  significant  problem. 

Using  a  100  ohm  standard  resistance  in  the  place  of 

the  equivalent  circuit  noise  was  measured  at  less 

-8 

than  10  volts  p.p.  on  the  X-Y  recorder,  for  a  circuit 

time  constant  of  0.3  seconds.  The  time  constant  is 

related  to  the  bandwidth  of  the  lock-in  output  by  the 

expression  B  =  1  . 

2t 
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CHAPTER  V 

EXPERIMENTAL  METHOD 

V. 1  Specimen  Preparation  and  Mounting  Procedure 

A  tunnel  junction  consists  of  a  metal-insulator- 
metal  sandwich  deposited  on  a  suitable  support  such  as  a 
glass  microscope  slide.  Junction  fabrication  was  begun 
by  thoroughly  cleansing  a  7/8"  x  3/8"  piece  of  glass 
plate  and  then  flaming  it  over  a  bunsen  burner  until  the 
edges  began  to  soften.  Indium  metal  contacts  were  then 
placed  on  the  glass  using  a  soldering  iron.  Leads  we re 
attached  to  the  contacts  and  then  the  substrate  was 
placed  above  a  mask  in  the  vacuum  deposition  bell  jar. 

Aluminum  was  used  to  form  the  normal  metal  in  this 
work.  This  was  deposited  on  the  slide  first.  A  vacuum 
of  ^2x10”^  mm  of  Hg  was  attained  before  any  evaporation 
was  begun. 

The  insulating  layer  in  the  junctions  used  was 
always  aluminum  oxide.  This  layer  was  formed  by  re¬ 
moving  the  substrate  from  the  vacuum  jar  after  the 
aluminum  evaporation  and  placing  it  into  an  oven  heated 
to  100°C  and  filled  with  atmospheric  air.  An  oxidation 
time  of  1  minute  usually  resulted  in  forming  enough  oxide 
on  the  aluminum  film  for  junction  resistances  of  between 
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100  and  300  ohms  after  the  Pb  strip  was  deposited  across 
the  oxide . 

Junction  areas  were  about  1 . 5  mm  square  and  their 
capacitance  was  normally  about  0.024  pf arads .  Two 
junctions  were  formed  on  each  substrate  so  that  six 
wires  were  required  for  four  terminal  measurement  of 
the  separate  junction  conductances.  This  doubled  the 
chance  for  a  successful  high  pressure  run. 

Usually  about  ten  attempts  were  required  before 
both  junctions  would  be  useable  on  the  same  substrate. 

The  specimen  resistance  and  capacitance  were 
measured  just  after  fabrication  using  a  bridge  built 
for  this  purpose  by  Rogers  (1964) .  They  were  then 
mounted  as  soon  as  possible  on  the  termination  strip 
screwed  to  the  bomb  cap.  Once  the  junctions  were  mount¬ 
ed  and  continuity  was  checked,,  the  bomb  chamber  was 
screwed  onto  the  nut  to  a  tightening  torque  of  about  300 
lbs -ft.  forming  a  seal  at  the  lense  ring.  The  experi¬ 
mental  chamber  inside  the  bomb  could  then  be  evacuated 
through  the  high  pressure  capillary.,  by  attaching  a 
vacuum  roughing  pump  at  an  outside  point.  Storage  in 
this  vacuum  greatly  increased  the  chance  of  retaining 
the  initial  characteristics  of  the  junction  while  the 
rest  of  the  steps  were  being  taken  to  commence  a  run. 

The  needle  valve  and  pumping  lines  were  then  attach- 
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ed  to  the  helium  reservoir.  The  thermometer  lead  con¬ 
nections  were  then  made  and  tested.  The  cryostat  vacuum 
can  was  then  soldered  into  place  using  Wood’s  metal 

alloy.  The  evacuation  of  the  cryostat  was  then  begun. 

_  q 

Once  the  vacuum  reached  ^2x10  mm  Hg*  the  roughing 
pump  was  removed  from  the  pressure  capillary.  The 
capillary  was  then  connected  to  the  high  pressure 
intensifier  section  and  the  whole  pressure  system*  in¬ 
cluding  the  working  chamber*  was  flushed  with  helium 
gas.  Once  this  was  done*  the  cryostat  was  cooled  to 
liquid  nitrogen  temperature.  Junctions  could  be  stored 
at  liquid  nitrogen  temperature  for  several  weeks  with 
constant  properties. 

V.2  Description  of  a  Typical  Run 

After  a  specimen  had  been  successfully  mounted  and 
tested*  evacuation  of  the  cryostat  was  begun  and  continued 

_7 

until  the  vacuum  was  better  than  5*10  mm  Hg  at  liquid 

_7 

nitrogen  temperature.  Often  it  was  better  than  1x10 
mm  Hg.  This  degree  of  evacuation  was  essential  because 
once  the  system  was  cooled  down  to  working  temperatures* 
the  residual  gas  would  constitute  a  heat  leak.  This 
might  then  result  in  a  temperature  gradient  across  the 
bomb . 

Before  commencing  the  cool-down  procedure  to  liquid 
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helium  temperature,  the  pressure  system  was  always  tested 
to  see  if  it  had  any  leaks.  This  test  was  usually  done 
to  pressures  of  about  55,-000  psi.  If  the  electrical 
lead  seal  held  and  none  of  the  cone  seals  leaked.,  the 
final  cool-down  to  liquid  helium  temperature  was  begun. 

Pressure  transmission  to  the  junctions  could  be 
checked  by  monitoring  the  junction  resistance  R^  as  a 
function  of  pressure.  This  resistance  decreased  sens¬ 
itively  as  the  oxide  is  compressed,  and  was  effectively 
used  as  a  means  of  measuring  pressure  at  the  junction. 
Capillary  blockage  during  pressurization  was  readily 
seen  since  suddenly  remained  constant  while  the  out¬ 
side  pressure  was  still  increasing.  The  changes  in  the 
junction  resistance  were  always  found  to  be  reversible 
with  pressure  change. 

V. 3  Low  Temperature  Pressure  Production 

The  helium  cool-down  procedure  was  started  by  empty¬ 
ing  the  inner  dewar  filled  with  liquid  nitrogen  and  re¬ 
installing  it  over  the  cryostat.  This  dewar  was  then 
sealed  to  the  recovery  system  so  that  the  He  gas  could 
be  collected  and  re-liquefied. 

A  larger  dewar  was  then  hung  around  the  inner  dewar 
and  filled  with  liquid  nitrogen.  Liquid  helium  from  a 
portable  storage  dewar  was  then  siphoned  under  a  He  gas 
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positive  pressure  head  into  the  inner  dewar .  This  process 

was  continued  until  the  system  was  cooled  down  to  about 
o  o 

30  -  35  Kelvin.  At  this  time  the  pressure  was  re¬ 

pumped  to  the  desired  working  value.  Heaters  on  the 
high  pressure  capillary  were  turned  on  to  keep  it  open. 
Cool-down  and  freezing  at  constant  pressure  was  then 
commenced  by  opening  the  needle  valve  and  sucking  helium 
liquid  from  below  the  vacuum  can  into  the  pumping  can 
mounted  on  the  pressure  bomb. 

Constant  pressure  freezing  of  56*000  psi  at  35°K 
resulted  in  ^50*000  psi  at  4.2°K.  The  pressure  vs 
curve  obtained  at  35°K  served  to  indicate  the  final 
frozen  in  pressure  at  the  working  temperature.  The 
temperature  variation  of  R^.*  for  zero  pressure*  is 
negligible  between  these  limits  so  no  correction  was 
made  for  the  pressurized  run  either. 

V. 4  Pressure  Dependence  of  the  Barrier 

The  tunneling  transition  probability  for  barrier 
penetration  as  used  in  equation  11-27  was  assumed  to  be 
independent  of  temperature  and  energy*  for  the  small 
energies  involved.  The  probability  of  barrier  penetra¬ 
tion  takes  the  form 

l  1  1 

I  Ml  =  exp-2/h/[2m(0(x)  -  E^)2  ]dx 

X2 


V-l 


% 

* 
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where  $(x)  is  the  potential  barrier  height  as  a  function 
of  x,  E  is  the  x-component  of  the  electron  kinetic 

energy,  and  m  is  the  electron  mass. 

2  1 

Assuming  that  |  M  |  oc  exp- (x^  -  x^))  where  pf  is 
an  average  potential  height,  and  x,  -  x^  is  the  barrier 
thickness,  pressurization  of  the  junction  should  increase 
P  approximately  exponentially  as  the  barrier  thickness 
decreases . 

Figures  11  and  12  show  resistance  vs  pressure  plots 
for  four  different  junctions.  Figure  13  shows  them 
plotted  on  a  log-linear  graph.  Notice  that  the  curves 
have  approximately  an  exponential  pressure  dependence  as 
one  would  expect.  The  lower  resistance  junctions  exhibit 
a  somewhat  larger  relative  change  for  a  given  pressure 
difference.  This  also  conforms  to  the  expectations  of 
an  idealized  theory. 

Deviations  from  a  true  exponential  behaviour  are 
probably  associated  with  effects  at  the  surfaces  sep¬ 
arating  the  oxide  from  the  metals. 

In  principle,  the  junction  capacitance  as  a  function 
of  pressure  should  also  be  obtainable.  In  practice, 
these  changes  were  not  large  enough  to  provide  much 
pressure  sensitivity. 

In  conclusion,  the  normal  junction  conductance 
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Figure  11 


Tunnel  resistance  vs  pressure 
for  junctions  1345B  and  1345. 
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Figure  12 


Tunnel  resistance  vs  pressure 
for  junctions  13^50  and  1345E. 
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Figure  13 


Semi- logarithmic  plot  of  tunnel 
resistance  vs  pressure. 
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is  sufficiently  pressure  sensitive  that  it  can  be 
used  as  a  pressure  monitoring  device  wherever  pressures 
of  a  few  thousand  bars  must  be  indirectly  measured  to 
moderate  accuracy. 
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CHAPTER  VI 

RESULTS  AND  DISCUSSIONS 

VI . 1  The  Pressure  Dependence  of  the  Transition  Temperature 
One  of  the  simplest  and  yet  one  of  the  severest  tests 
of  the  pressure  method  used  in  this  work  concerns  the 
pressure  dependence  of  the  critical  temperature  of  Pb . 

This  is,  therefore.,  discussed  first. 

The  Pb  transition  temperature  was  obtained  for  each 
junction  by  monitoring  the  zero  bias  conductance  while 
the  temperature  was  slowly  decreased.  At  the  transition 
temperature,  the  energy  gap  began  to  develop  and  this 
appeared  as  an  abrupt  decrease  in  a.  The  transition  was 
well  defined  in  all  junctions  measured  and  an  average 
value  of  four  readings  constitutes  each  quoted  tempera¬ 
ture.  While  an  error  estimate  of  +  7  mdeg  is  placed  on 
the  absolute  temperature  near  7°K,  an  estimate  of  +  4  mdeg 
error  is  placed  on  the  difference  AT  due  to  the  appli- 
cation  of  a  pressure  (measured  in  bars) .  Table  1  shows 
the  shift  in  transition  temperature  for  several  junctions. 

While  the  bulk  Pb  transition  temperature  is 
7d93  +  0.005°K,  (Franck  and  Martin,  1961),  most  of  the 
thin-film  zero-pressure  Pb  transition  temperatures  were 
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slightly  higher  than  this  value.  This  effect  is  assoc¬ 
iated  with  thin  film  geometry,  and  is  present  in  most 
superconductors  deposited  as  films. 


TABLE  1 


JCN 

Tc(P=0) 

T0(P) 

AP 

AT 

c 

• dTc/dP 

°K 

°K 

bar 

°K 

°K/b; 

ar 

1345A 

6.965 

6.885 

2378o0 

0.080 

3.36 

+ 

0.33 

X 

10'5 

1345B 

7.246 

7.158 

2551 o2 

0.088 

3.45 

+ 

0.32 

X 

10~5 

1246b 

7.240 

7.197 

1309.9 

0.053 

4.05 

+ 

0.50 

X 

10"5 

1345C 

7.240 

7.110 

3447.5 

0.119 

3.45 

+ 

0.32 

X 

10"5 

1246c 

7.248 

7.140 

3172.0 

0.108 

3.40 

-f 

0.32 

X 

10"5 

1345D 

7.225 

7.125 

2757.2 

0.100 

3.63 

+ 

.3- 

00 

• 

0 

X 

10-5 

1246d 

7.225 

7.120 

2757.2 

0.105 

1 - 1 

CO 

• 

CO 

+ 

0.35 

X 

10"5 

1345E 

7.243 

7.142 

2730.0 

0.101 

3.69 

+ 

0.35 

X 

10-5 

Average  = 

3.61 

t 

T 

0.35 

X 

10"5 

A  pressure  uncertainty 

'  of  ^6/0 

is  included  in 

the  error 

estimate . 

Table  2  shows  a  list  of  published  values  of  dTc/dP 
and  dH  /dP  for  Pb .  All  of  the  values  of  dT  /dP  (except 
the  value  measured  during  this  work)  were  obtained  from 
magnetic  field  measurements  on  bulk  samples.  This  can 
be  done  by  utilizing  the  relation, 

dT  oh /ap) 

- -  - - - - =•  ,  T-T  .  VI-1 

dP  ( aHc/c>T)  p 


s.  * 

> 

' 
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The  agreement  between  thin  film  and  bulk  values  is 
quite  close  and  generates  confidence  in  the  pressure 
technique.  Obviously,  if  strain  effects  are  present  due 
to  differential  contraction  between  the  glass  substrate 
and  the  films,  their  effects  are  not  large  enough  to  pose 
a  major  problem. 


TABLE  2 


(SH  /SP)T 

C 


( dTc/dP )  Ez^0  AUTHORS 


(Gauss/bar) 


( °K/bar ) 


-9.62  +  0ol8  x  10  3 


-9.03  +  0.15  x  10  3 


-4.04  +  0.08  x  10"5 
-3.86  +  0.12  x  10'5 
-3.79  +  0.07  x  10~5 

-3 060  +  0.35  x  10"5 


R.R.  Hake  and 
D.E.  Mapother 
(1956) 

T.F.  Smith  and 

C. W.  Chu 
(1967) 

M.  Garfinkel  and 

D. E.  Mapother 

(1961) 

This  work 
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VI . 2  Phonon  Shifts  in  Pb  Under  Pressure 

In  section  II. 8  a  discussion  was  given  showing  how 
the  predominant  phonon-peaks  are  reflected  as  structure 
in  the  density  of  states  (a)  for  strong-coupling  Pb . 

The  phonon  emission  process  is  the  mechanism  leading  to 
this  structure ^  and  the  structure  is  of  such  a  form  that  a 
passes  through  an  inflection  point,,  while  da/dV  passes 
through  a  minimum  at  the  predominant  phonon  energies. 

These  energies  occur  near  the  ends  of  the  transverse  and 
longitudinal  phonon  branches. 

We  have  measured  the  shifts  in  a  and  da/dV  at 
and  cd £  for  ten  junctions  subjected  to  hydrostatic  pressure. 
Figure  14  shows  conductance  results  taken  at  1.4°K  for 
a  typical  junction.  Since  the  X-Y  plotter  paper  is 
14"  x  17" ,  the  plot  has  been  photographically  reduced  for 
presentation  here.  The  phonon  energies  are  measured  from 
the  gap-edge  energy ^  so  the  P  =  0  and  P  >  0  traces  have 
been  re-drawTn  onto  a  single  graph  with  the  gap-edge 
energies  set  to  coincide.  Figures  15  and  16  show  actual 
size  traces  of  a  and  da/dV  for  the  transverse  and  long¬ 
itudinal  regions  respectively.  These  figures  display 
several  results: 

(1)  The  shift  in  2Aq  is  directly  observable  at  the  left 
in  figure  14. 

(2)  The  shift  in  is  easily  resolved  in  both  a  and  da/dV. 
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Figure  14 


The  normalized  conductance 
1345C  at  P  =  0  bar  and  P  = 
gap-edge  energies  are  set 


a  for  junction 
3172.0  bar.  The 
to  coincide. 
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Figure  .15 


Normalized  conductance  a  and  second  derivative 
2  2 

d  i/dv  (in  arbitrary  units)  for  junction  1345C 
at  T  =  1.4°K.  The  energy  range  near  the  longi¬ 
tudinal  phonon  peak  is  shown.  Solid  lines  are 
at  P  =  Oj  broken  lines  are  at  P  =  3445  bar. 
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Figure  1 6 


Normalized  conductance  a  and  second  derivative 
2  2 

di/d v  (in  arbitrary  units)  for  junction  134pC 
at  T  =  1.4°K.  The  energy  range  near  the  trans¬ 
verse  phonon  peak  is  shown.  Solid  lines  are  at 
P  =  0,  broken  lines  are  at  P  =  3445  bar. 
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(3)  There  is  a  general  structure  amplitude  reduction 

p 

with  pressure  associated  with  changes  in  a^(ai)F(a))  and 
the  density  of  states „ 

(4)  Because  of  (3)*  the  shift  in  co^  is  not  as  clearly 

defined  as  for  the  longitudinal  case.  Also,  because  of 

the  separate  contributions  of  cd,  ,  and  cjd.  ,  dc/dV  has 

_  1  2 

some  structure  at  which  renders  its  use  impossible  for 
measuring  pressure  shifts.  A  new  criterion  for  measuring 
oxj.  was  therefore  established.  Its  value  for  a  particular 
curve  was  obtained  from  the  mid-height  (inflection)  point 
in  the  conductance  structure  drop-off.  Table  3  lists  a 
number  of  values  of  Scd^/SP  obtained  in  this  way.  Table  4 
lists  values  of  Sa^/SPo 

TABLE  3 


JUNCTION 

o>t  ( P=0 ) 

o>t(P>0) 

SP 

Scot/6P 

meV 

meV 

bar 

meV/bar 

1345A 

4 . 585 

4  0630 

1930.6 

2.33  +  1.00 

X 

10  5 

13^55 

4.580 

4.650 

2551.2 

2.74  +  1.00 

X 

i— 1 

0 

1 

\j 1 

1345C1 

4.575 

4.642 

3447.5 

1.94  +  1.00 

X 

10"5 

I345C2 

4.575 

4.650 

3447.5 

2.18  +  0.60 

X 

10"5 

1345C3 

4.585 

4  „  660 

3172.0 

2.36  +  0.50 

X 

10'5 

1345D 

4.600 

4.675 

2757.2 

2.72  +  1.00 

X 

10"5 

The  average  for  co^(P=0)  =  4.585  meV,  and  the  average  for 
6mt/6P  is  2.37  +  °-7°  x  10’5 


meV/bar.  Therefore  we  have. 


, 

* 
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dikicD 

=  5.17  + 

1.50  x  10 

-6  ,  -1 

bar 

dP 

TABLE 

4 

JUNCTION 

^(p=o) 

u/P>0) 

SP 

6cd£/6P 

meV 

meV 

bar 

meV/bar 

13^5A 

8.710 

8.825 

1930.6 

5.96  +  1.20 

X 

10“5 

1345B 

8.650 

8.810 

2551.2 

6.27  +  1.20 

X 

10"5 

1345C1 

8.650 

8.835 

3447.5 

5.37  +  0.80 

X 

10"5 

1345C2 

8.675 

8.863 

3447.5 

5.45  +  0.80 

X 

10'5 

1345C3 

8.725 

8.930 

3172.0 

6.46  +  1.00 

X 

10'5 

1345D 

8.700 

8.863 

2757.2 

5.91  +  1.00 

X 

H 

O 

l 

vn 

average  (P=0)  - 

8.685  me 

V,  and  the 

:  average  for 

6a>£/6 

Is  5.90  +  1.00  x  10  J  meV/bar.  Thus,  In  this  case  we 
have 


dinco 


dP 


1  =  6.80  +  1.00  x  10  6  bar  1 


Since  the  minima  in  da/dV  occur  at  peaks  in  the  phonon 
spectrum,  a  shift  in  dc/dV  as  the  superconductor  is  sub¬ 
jected  to  hydrostatic  pressure,  gives  directly  the 
Gruneisen  y  for  the  predominant  modes  through  the  expres¬ 
sion 

dina) .  6cd  ./cd  . 

_ J  =  _ B 


VI-2 


dinV 


6P 


where  B  is  the  bulk  modulus,  and  6co .  is  the  shift  in  mode 
j  due  to  a  pressure  P.  Using  the  value  of  B  =  4.88  x  10^ 


■ 


82 


bar,  due  to  Waldorf  and  Alers  (1962),  gives  the  Gruneisen 
values 


=  (6.80  +  1.00) (0.488)  =  3.32  +  0.49 
7t  =  (5.17  +  lo50) (0.488)  =  2.52  +  0.74 


This  result  is  true  for  the  predominant  modes,  despite 
a  possible  change  in  the  electron-phonon  coupling  constant 
(see  Appendix  A). 

It  is  instructive  at  this  stage  to  briefly  outline 
the  general  Gruneisen  7  theory  in  order  to  facilitate 
comparison  of  our  results  with  results  obtained  by  other 
workers.  The  Gruneisen  rule  states  that  at  low  temper¬ 
ature,  the  coefficient  of  thermal  expansion  a  is  propor¬ 
tional  to  the  constant  volume  specific  heat,  C  ,  through 
the  relation, 

Va 


VI-3 


where  K  is  the  compressibility.  This  result  follows 
strictly  when  all  7.  given  by 

J 

dinv  . 

7  - - 1  VI-4 

J  dinV 


t  h 

are  equal.  The  subscript  j  refers  to  the  j  normal 


mode  of  vibration,  and  is  related  to  these  mode  gammas 


by 


7G  =  J 


2  7  .C  . 

vj  . 


2  C. 


vj 


VI-5 


. 
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C  .  is  the  Einstein  heat  capacity  of  the  mode. 

v  J 

Blackman  (1957*  1959)  has  derived  a  second  form 
under  which  VI-3  is  also  strictly  valid.  This  depends 
on  forming  a  quantity  y( v)  which  is  the  average  of  y(v) 
taken  between  neighbouring  frequency  contours.  Then., 


/  7(v)p(v)Cvvdv 


= 


J  p(v)Cvvdv 


VI-6 


holds  strictly  when  y(v)  =  constant.  Here  p(v)  is  the 
frequency  distribution  function.  Blackman  has  shown 
that  this  is  true  for  two  special  cases: 

(1)  At  sufficiently  high  temperatures,  the  heat  capacity 
of  each  mode  is  constant  =  3/2  k. 

o 

(2)  In  the  region  where  the  law  for  Cv  and  cl  holds 
sufficiently  well,  can  be  derived  from  the  free  energy 
and  takes  the  form 


III  Z  c’3y^  d9 
,1=1  J  J 

3  2 

HI  2  c'3  d9 

J=1  J 


VI-7 


where 


7  j  =  1/3  -  V/c^dOj/dV) 


vi-8 


and  c  .  is  the  velocity  of  a  wave  in  a  given  direction 
J 

in  one  of  the  three  acoustic  branches  while  V  is  the 
volume.  > For  cubic  crystals  it  can  be  shown  that  the 
velocity  of  the  [100]  longitudinal  plane  wave  is  (  cu/p> 


. 

ij  ,%q  *  iev  i.  Si’  ri.'i  r.4  ,  &  &■  ( 
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where  p  is  the  density.  At  the  same  time,  the  velocity 


of  the  transverse  wave  is  given  by  (C///p)2.  In  these 

expressions,  C-^  and  are  lattice  constants  associated 

with  propagation  in  the  continuum  region  of  the  dispersion 

0 

curve.  Thus,  using  the  definition  of  y.  one  has. 


y^[ 100]  =  1/3  -  din(cn/p)2 

dinV 

=  1/3  -  1/2 


dinC-,-,  ,  dfnp 


11 


1/2 


VI-9 


VI-10 


dinV 


dinV 


Since  =  -  ^  =  -1, 

dinV  dinV 

7^100]  =  -  1/2  dinCll  -  1/6.  VI-11 

dinV 

Similarly, 

7t[l00]  =  -  1/2  dinC44  -  1/6.  VI-12 

dinV 


Acoustic  measurements  yield  the  velocity,  and  hence  the 
elastic  constants  and  y  for  a  given  direction  of  sound 
propagation  in  the  non-dispersive  region.  For  a  general 
comparison  with  other  fee  metals,  y’s  measured  by  acoustic 
experiments  (i.e.  continuum  results)  are  given  in  Table  5 
along  with  the  present  result. 

There  seems  to  be  very  little  experimental  data 
available  on  Pb  for  comparison  with  this  work.  Phonon 
dispersion  curves  have  been  measured  in  Pb  at  100°K  by 
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TABLE  5 


DIRECTION 

0 

°D 

7f 

7t 

rl 

7t 

Z2 

At .No. 

AUTHOR 

Na 

[100] 

150 

1.51 

1.06 

1.06 

11 

W .B. Daniels 

[no] 

1.36 

1.06 

1.06 

(i960) 

[in] 

1.34 

1.06 

1.06 

Al 

[100] 

390 

2.28 

2.80 

2.80 

13 

R.E.Schmunck  & 

[no] 

2.43 

2.80 

2.36 

Charles  S . Smith 

[mi 

2o43 

2.53 

2.53 

(1959) 

Cu 

[100] 

315 

2.48 

1.92 

1.92 

29 

W.B. Daniels  & 

[no] 

2.30 

1.92 

1.49 

Charles  S. Smith 

[in] 

2.19 

1.76 

1.76 

(1958) 

Ag 

[100] 

215 

2.71 

2.38 

2.38 

47 

W.B. Daniels  & 

[no] 

2.69 

2.38 

1.96 

Charles  S. Smith 

[in] 

2 . 68 

2.21 

2.21 

(1958) 

Au 

[  100  ] 

170 

2.86 

3.38 

3.38 

79 

W .B .Daniels  & 

[no] 

3.00 

3.38 

2.31 

Charles  S. Smith 

[in] 

3.03 

2.94 

2.94 

(1958) 

Pb 

Av . 

105 

3.32 

2.52 

2.52 

82 

This  work 

Brockhouse  et  al  (196IA)  using  inelastic  neutron  scatter¬ 
ing.  They  have  also  measured  the  frequency  shift  due  to 
change  in  temperature  from  100°K  to  425°K  for  a  few  phonon 
modes  (196IB).  More  recently,  Stedman  et  al  (1967A),  again 
using  inelastic  neutron  scattering,  have  made  measurements 
at  80°K  and  300°K.  The  total  phonon  spectrum  F(cjd)  vs  gd 
requires  the  summation  of  the  dispersion  curves  over  all 


• 
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directions  and  is,  therefore.,  not  easily  obtainable  c  At 

present,  F(ai)  is  not  known  as  a  function  of  temperature. 

However ,  Stedman  et  al  have  succeeded  in  calculating  F(a>) 

vs  a)  for  Pb  at  80°K  (1967B).  They  used  their  dispersion 

results  for  three  measured  symmetry  directions  of  the 

three  main  branches.  Values  for  intermediate  directions 

were  computer  interpolated,  and  therefore  the  resultant 

spectrum  is  obtained  without  making  any  assumptions 

based  on  models  of  the  lattice.  Their  phonon  spectrum 

is  almost  identical  to  the  published  results  of  McMillan 

o 

and  Rowell  (figure  2)  on  a  (a))F(a>)  vs  cd  obtained  by 
electron  tunneling 0 

Gilat  (1965)  has  also  calculated  the  Pb  phonon 
spectrum  from  a  general  Born-Von-Karman  model  and  know¬ 
ledge  of  the  dispersion  curves  of  Brockhouse  et  al .  His 
result  is  expected  to  be  less  accurate  since  it  is  pri¬ 
marily  a  mathematical  model  calculation,  and.,  in  fact,  it 
agrees  poorly  with  the  above  mentioned  results. 

While  the  neutron  results  are  interesting  in  that 
they  show  generally  good  agreement  with  the  tunneling 
phonon  spectrum,  no  information  regarding  the  Gruneisen 
y’s  is  given,  and  these  are  presently  under  investigation. 
Lechner  and  Quittner  (1966)  have  measured  pressure- 
induced  phonon  frequency  shifts  for  six  true  phonon  modes 


. 

, 
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in  Pb  at  room  temperature.,  and  3000  atmospheres  pressure. 
They  used  inelastic  neutron  scattering  for  these  measure¬ 
ments,  and  the  G-runeisen  y!s  for  these  modes  are  listed 
in  Table  6. 


TABLE  6 


DIRECTION 

POLARIZATION 

aq/27T 

7  • 

[  Hi] 

T 

0.867 

2.66  + 

0.30 

[ 

L 

0.867 

0.91  + 

0.43 

Uoo] 

T 

1.000 

3.86  + 

0.78 

Uoo] 

T 

0.650 

0.97  + 

0.62 

Uoo] 

L 

1.000 

1.07  + 

1.39 

Uoo] 

L 

0.750 

2.29  + 

1.31 

Although  Townsend  (1964)  has  indicated  that  Pb  crystal¬ 
lites  deposit  preferentially  with  their  (ill)  axes  aligned 
perpendicularly  to  the  film  surface.,  the  attitude  of  Rochlin 
(1967)*  that  tunnel  results  yield  a  directionally  averaged 
phonon  spectrum.,  is  taken  by  this  author.  Therefore  shifts 
in  cd^_  and  cd^  obtained  in  this  work  cannot  be  directly  com¬ 
pared  to  the  neutron  diffraction  results  of  Lechner  and 
Quittner . 
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Payne  (1964)  was  the  first  to  publish  results  in 
which  the  tunneling  technique  was  used  to  study  shifts 
in  phonon  modes  under  an  applied  stress.  He  measured 
shifts  in  some  of  the  predominant  phonon  modes  in  german¬ 
ium  subjected  to  uniaxial  compression,  but  he  notes  that 
considerable  care  must  be  taken  in  accounting  for  shear 
induced  effects  resulting  from  this  form  of  "pressure" 
production.  This  problem  is  absent  when  the  junctions 
are  pressurized  in  a  solid  helium  bath. 

Hodder  and  Briscoe  (1966)  have  attempted  to  measure 
the  shifts  in  the  phonon  structure  in  Pb,  by  bending 
junctions  of  Pb-Pb-oxide-Pb,  at  liquid  helium  temperatures. 
These  junctions  exhibit  structure  in  the  neighbourhood 
of  the  transverse  phonon  peak  which  can  be  identified  with 
the  Van  Hove  singularities  in  the  phonon  spectrum.  The 
nature  of  these  singularities  has  been  discussed  by  Van 
Hove  (1953).,  by  Rosenstock  (1955).,  and  by  Phillips  (1956), 
while  their  presence  in  tunneling  results  has  been  explain¬ 
ed  by  Rowell,  Anderson,  and  Thomas  (1963),  and  by  Scalapino 
and  Anderson  (1964) .  Hodder  and  Briscoe  therefore  shift 
these  singularities  by  bending  their  junctions.  However, 
since  there  are  several  singularities  in  the  energy  region 
of  od^_ ,  and  since  they  are  not  completely  resolved  from 
each  other,  it  is  doubtful  if  their  measured  shifts  corres- 
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pond  to  actual  phonon  shifts.  Also,  their  results  depend 
crucially  on  the  choice  of  Poisson's  ratio,  since  some 

effort  was  made  to  calculate  shear  effects.  This  can 

best  be  illustrated  by  their  calculated  value  of  6fnV/unit 

strain,  which  varies  from  0.229  -  0.519*  depending  on 

the  assumptions  me.de.  On  top  of  these  uncertainties  is 

another  large  contribution  due  to  the  extremely  small 

phonon  shifts  produced.  They  find,  in  this  way,  that 

6fnar /unit  strain  ^  2  6incD  ,/unit  strain.  This  is  in 

disagreement  with  our  results. 

The  thermal  expansion  measurements  due  to  White 

(1962),  possibly  yield  the  best  comparison  for  the  present 

results.  If  one  assumes  that  y(ax  )  and  7(05.)  remain 

r  jo 

constant  over  the  whole  energy  range  of  the  two  branches, 
then  at  low  temperature,  where  the  transverse  modes  are 
predominantly  excited,  equation  VI-7  becomes 

2/3(7tA  l)  +  1/3(7 /<^) 

7(J  2/3(1/ C3  )  +  1/3(1/ c^) 

7t  +  l/2(7i) (cfc/c^)3 
1  +  l/2(ct/Ci)3 

~  7^-  =  2.52  +  0 .73 . 

At  high  temperature,  where  all  the  modes  are  excited. 


90 


one  has, 

7q  =  2/3  7(o>t)  +  1/3  7(a>i) 

=  2.78  +  0.64. 

These  values  compare  with  by  White  of  2.6  and  2.7 
respectively  (no  uncertainty  limits  are  given) .  While 
the  present  results  seem  to  show  more  variation  than  the 
thermal  expansion  results*  the  two  sets  of  values  probably 
agree  within  experimental  error. 


.  ’  •  ,  ' 
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VI . 3  The  Pressure  Dependence  of  the  Energy  Gap 

Energy  gap  values  were  obtained  as  a  function  of 
temperature  by  relating  the  zero  bias  conductance  to 
tabulated  values  of  Aq  found  in  Bermon 1 s  (1964)  table. 
His  calculations  relate  the  normalized  density  of 
states  (a)  to  Aq  for  a  BGS  constant  gap-parameter 
theory.  Although,  strictly  speaking,  the  gap  at  the 
gap-edge  energy  should  have  been  measured,  the  use  of 
the  conductance  minimum  gives  a  more  consistent  method 
of  estimating  A  (T),  especially  when  making  the  pressure 
comparison.  Furthermore,  since  A^(co)  ~  constant  and 
A^cd)  rv,  0  for  energies  up  to  the  gap-edge  energy  and 
a  little  beyond,  one  might  expect  the  temperature  vari¬ 
ation  of  the  gap  to  follow  a  BCS  behaviour  to  a  good 
approximation.  We  therefore  assume  that  shifts  in  the 
zero  bias  gap  value  are  in  1:1  correspondence  with 
shifts  at  the  gap-edge  energy. 

Figure  17  shows  a  plot  of  A  (T)  vs  T  for  Junction 
1345E.  The  lower  curve  was  taken  for  P  =  39*600  psi 
while  the  upper  consists  of  P  =  0  points  taken  before 
and  after  the  pressure  run.  The  reproducibility  of  the 
P  =  0  points  rules  out  permanent  sample  deformation  as 
a  possible  origin  of  the  effect.  Points  below  4.2°K 
become  increasingly  less  accurate  because  a  is  now 
close  to  zero  and  is  only  slightly  temperature  dependent. 


Figure  17 


Energy  gap  (at  the  gap  edge)  of  Pb  as  a  function 
of  temperature.  Order  of  runs:  black  dot  (P  =  0) 
open  circles  (P  =  2730  bar);  black  squares  (P  =  0 


. 
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Figure  18  shows  a  plot  of  (A  (T)/kT  vs  t  for 

O  0 

the  same  junction.  The  BCS  weak- coupling  limit  predic¬ 
tion,  taken  from  calculations  by  Muhlschlegel  (1959)* 
is  plotted  as  well.  The  fact  that  the  P  =  0  and  P  >  0 
plots  do  not  coincide.,  is  direct  evidence  that  the 
ratio  is  pressure  dependent. 

Figure  19  shows  the  same  points  plotted  in  the  form 
(A  /kT  )/(A  /kT  )tdPQ  vs  t  for  the  two  pressures.  Both 

O  C  O  C  -DUkD 

sets  of  points  lie  on  essentially  straight  lines  through 
10266  and  102480  Although  the  points  for  t  >  0.8  seem 
to  lie  somewhat  lower  than  the  values  for  t  <  0.8,  be¬ 
cause  of  the  experimental  uncertainty,  it  was  decided 
that  this  deviation  should  not  be  considered  significant. 
One  may  conclude  from  this  plot,  in  agreement  with  the 
predictions  of  Scalapino,  Swihart,  and  Wada  (section  II. 9)* 
that  the  gap  vs  temperature  relation  follows  a  BCS  type 
variation  to  "within  a  few  percent" .  This  seems  to  be 
true,  even  under  pressure. 

From  the  difference  in  the  least  squares  fitted 
horizontal  straight  lines  through  the  P  =  0  and  P  >  0 
plots,  one  obtains, 

2A 

6(— )t=const  =  (1.248  -  1.266)  (3.528)  =-0.0635. 

kT 


i- 
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Figure  18 


Squared  reduced  energy  gap  (at  the  gap  edge)  of  Pb 
as  a  function  of  reduced  temperature.  Black  dots 
and  squares,  P  =  0;  open  circles,  P  =  2730  bar. 


t-T/Tc 


Figure  19 


Plot  of  ( A  /kT  )/( A  /kT  vs  reduced 

temperature;  black  dots  and  squares,  P 
open  circles,  P  =  2730  bar. 
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1.35 


o 


i/i  = 


96 


Thus, 

din(2Ao/kTc) 

dP 


0.0635/4.403 

-  =  -(5.3  +  0.5)  x  10  oar 

2730 o 5 


and  from  the  relation 


din(2A  /kT  )  dinA  dinT 
v  o'  c '  o  c 


dP 


dP 


dP 


dinT 


we  have,  using 


dP 


c  =  -5.06  +  0 .30  x  10  6bar  1 


VI-13 


d inA  _r  , 

- -=  -10.3  +  1.1  x  10  Dbar 

dP 

Thus  A  shifts  with  pressure  almost  twice  as  fast  as  T 
o  *  c 

does . 

Additional  Information  regarding  the  shift  In  Aq 
was  obtained  by  measuring  the  shift  In  the  2°K  gap  wall 
as  a  function  of  pressure.  At  a  temperature  of  2°K  the 
gap  Is  99.28/  fully  developed  so  that  we  may  assume 


dinA  (0)  din A  (2°K) 

o  _  _ o_ _ _ 

dP  dP 


VI-14 


Figures  20  and  21  show  typical  results  for  two  junctions. 
The  values  are  listed  In  table  7.  Since  the  shift  in  t 
at  2°K  due  to  the  applied  pressure  is  negligible,  one  can 


m 


. 
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Figure  20 


The  gap  region  conductance  at 
junction  13450.  Pressures  are 


and  P  =  2730  bar. 


mV 


Figure  21 


The  gap  region  conductance  at  2° 
junction  13450^ .  Pressures  are., 


C  fo 
P  = 


and  P  =  3^45  bar. 


. 


99 


TABLE  7 


JUNCTION 

1345C 

1345D 

1345E 

Average 


dinAo  bar"1 


dP 


-9.52  +  0.80  x  10 
-9.60  +  0.80  x  10 
-9.65  +  0.80  x  10 
-9.6  +  0.8  x  10 


dinT  .  Q  -1 
c  oar 


dP 


-4. 8l  +  0.50  x  10 

-4.79  +  0.50  x  10 

-5.02  +  0.50  x  10 

-4.9  +  0.5  x  10 


again  refer  to  dinA  /dP  as  occurring  at  constant  reduced 

o 

temperature.  Substituting  average  values  of  dinA  /dP  and 

dinT  /dP  from  table  7  yields: 
c 


dfn(2A  /kT  )  r 

-  =  -(4.7  +  1.3)  x  10  bar 

dP 


Because  of  the  experimental  uncertainty,  the  2°K  shift 
is  not  considered  to  be  significantly  different  from  the 
extended  temperature  shift  measured  for  Junction  1345E. 
Taking  an  average  value  to  represent  the  experimental 
measurements  gives. 


din ( 2A  /kT  )  _6 

- - - —  =  -(5.0  +  0.9)  x  10  bar  -1 

dP 


In  the  BCS  weak- coupling  theory  both  Tc  and  A  (0)  have 


' 
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the  form  Tc°c  Aqc c  0^exp(-l/N(O) V)  .  Thus  pressurizing  the 
system  results  in  identical  shifts  in  both  T  and  A  (0) 
so  that  the  ratio  2AQ(0)/kTc  should  remain  constant. 
However,  for  the  strong-coupling  ratio  measured  for  Pb, 
no  such  simple  relation  exists  between  the  two  parameters. 
We  guess  that  the  pressure  dependence  of  the  Pb  ratio  has 
its  origin  in  a  change  in  the  electron-phonon  interaction 
strength.,  along  with  changes  in  the  density  of  states. 
McMillan  (1968)  has  published  an  abstract  of  a  paper  to 
appear  in  Phys .  Rev.;  in  which  he  finds  the  coupling 
constant  is  strongly  frequency  dependent;  the  coupling 
strength  decreasing  as  the  phonon  frequencies  are  in¬ 
creased  0  Thus  by  changing  F(cd)  to  higher  values  with 
pressure;  the  electron-phonon  interaction  is  somewhat 
reduced  toward  a  BCS  limit  as  is  observed  in  the  gap/Tc 
ratio . 

This  argument  is  further  substantiated  by  the  work 
of  Zavaritskii  (1967).  He  finds  that  by  depositing  the 
Pb  for  his  Pb-Pb  junctions  at  1.6°K;  a  maximum  amount  of 
crystalline  lattice  strain  is  produced.  The  associated 
phonon  spectrum  exhibits  many  more  low  energy  modes;  and 
since  these  modes  contribute  most  in  the  gap  parameter; 

Aq  is  increased.  Subsequent  annealling  of  his  samples 
from  100°K  to  room  temperature  in  stages,  increases  the 

o 

definition  in  a  (a))F(a)),  and  the  spectrum  gradually 
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reverts  back  to  its  usual  form  with  an  accompanying  de¬ 
crease  in  the  gap-parameter.  He  does  not  measure  the 
change  in  the  transition  temperature  as  a  function  of 
annealling,  so  nothing  can  be  said  about  the  gap/Tc  ratio. 
However.,  it  appears  that  the  frequency  dependent  nature 
of  the  coupling  constant  is  further  borne  out . 

Although  this  author  feels  confident  that  this  ex¬ 
planation  of  the  measured  pressure  dependence  of  A  (0)/kT 

O  0 

is  correct.,  two  sources  of  possible  anomalous  behaviour 
cannot  be  completely  discounted.  These  include  aniso¬ 
tropy  and  low  temperature  strain  arguments. 

Anisotropy  of  the  energy  gap  in  Pb  has  been  studied 
by  Richards  and  Tinkham  (i960),  by  Townsend  and  Sutton 
(1963),  and  by  Rochlin  (1967)0  Richards  and  Tinkham  per¬ 
formed  electromagnetic  absorption  experiments,  while  the 
latter  two  investigations  were  carried  out  with  tunnel 
junctions.  In  the  tunnel  junctions,  the  anisotropy,  which 
leads  to  two  or  more  gap  values,  is  a  function  of  film 
thickness  and  disappears  for  films  less  than  ^0.7  H  in 
thickness.  In  these  thin  films,  due  to  the  small  electron 
mean  free  path,  the  superconductor  takes  on  the  behaviour 
of  a  "dirty"  superconductor  and  exhibits  a  single  average 
energy  gap,  even  for  high  purity  samples.  (The  Pb  and  A1 

films  used  in  the  present  work  for  junction  fabrication, 

o  o 

are  between  1,000  A  and  3^000  A  in  thickness,  so  aniso¬ 
tropy  effects  should  not  be  too  significant.) 
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Two  schools  of  thought  exist  concerning  the  origin 
of  the  gap  anisotropy.  Townsend  and  Sutton  ascribe  it  to 
Fermi  surface  irregularities,  while  Bennett  (1965)  and 
Rochlin  (1967)  feel  it  is  associated  with  anisotropy  in 
the  phonon  spectrum. 

Attempts  to  associate  the  pressure  dependence  of 
2A  (0)/kT  with  an  increase  in  gap  anisotropy  (due  to 
Fermi  surface  distortion),  can  be  ruled  out,  in  view  of 
the  work  done  on  the  pressure  dependence  of  the  Fermi 
surface  of  Pb  „  Anderson  and  Gold  (1965)*  and  Anderson 9 
O'Sullivan,  and  Schirber  (1967).,  have  carried  out  de  Haas- 
Van  Alphen  measurements  on  Pb  subjected  to  hydrostatic 
pressure o  In  the  first  paper,  the  Fermi  surface  is  found 
to  scale  in  all  directions  under  pressure  at  a  rate  close 
to  free-electron  model  predictions.  In  the  second  paper, 
this  is  somewhat  revised  in  that  isotropic  scaling  still 
occurs,  but  now  it  scales  at  approximately  2  times  free 
electron  predictions.  In  any  event,  excessive  distortion 
of  the  Fermi  surface  due  to  pressurizing,  seems  to  be  an 
unlikely  explanation  of  the  observed  pressure  dependence 
of  the  ratio  2A  (0)/kT  0  The  association  of  the  result 
with  increased  phonon  anisotropy  under  pressure  remains 
as  a  possible  explanation  of  the  effect  in  terms  of  gap 
anisotropy.  If  an  average  gap  value  is  measured,  with 
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directional  contributions  from  relations  of  the  form 
2Aq1(0)  -  m1  kTcl^  2Aq2(0)  =  mg  kTQ2,  .  .  .  etc,,  since 
only  the  maximum  transition  temperature  value  T  is 

measured,  the  result  dA  n  (0)/dT  4-  constant  implies 

that  either  the  rrn  values  are  pressure  dependent  (the 
explanation  favoured  by  this  author)  or  that  some  of 
dAQ2(0)/dP,  dAo2(0)/dP,  .  .  .  are  not  equal.  The  only 
sure  test  of  this  would  be  to  pressurize  Pb-Pb  junctions 
with  detectable  gap  anisotropy  and  try  to  relate  dT  /dP 
to  the  pressure  dependence  of  the  largest  gap  value. 

The  remaining  possibility  for  explaining  the  pressure 
dependence  of  2A  (0)/kT  is  film  strain  at  low  temperature. 
Caswell,,  Priest,  and  Eudo  (1963)*  have  measured  low  temp¬ 
erature  strain  in  Pb  films  and  find  that  the  amount  depends 
on  the  evaporation  conditions ,  the  amount  of  annealling 
carried  out  at  room  temperature *  and  on  film  thickness. 

This  therefore  parallels  (and  precedes)  the  results  obtain¬ 
ed  by  Zavaritskii  (1967).  All  of  the  junctions  used  in 
the  present  work  were  annealled  at  room  temperature  for 
at  least  one  day.,  and  sometimes  for  up  to  three  weeks. 

Since  the  results  are  all  similar.,  the  annealling  time 
does  not  seem  to  effect  the  pressure  dependent  result. 

This  is  probably  due  to  the  fact  that  differential  effects 
are  measured.  Also,  the  annealled  condition  of  a  given 
sample  can  be  considered  constant,  since  no  temperature 


■  ■ 


104 


annealling  occurs  until  the  temperature  exceeds  100°K 
(Zavarltskli  (1967)).  Between  the  P  =  0  and  P  >  0  runs 
the  temperature  never  goes  above  *\j  40°K. 

Wu  (1967)*  Adler,  Jackson,  and  Will  (1967)*  and 
Claeson  (1966,  1967),  have  changed  the  electron-phonon 
coupling  strength  in  Pb  by  alloying  with  indium  and 
thallium „  Their  results  have  been  summarized  in  Table  8 
(taken  from  Wu) . 

TABLE  8 


Alloy  Electrons 

Ao(0) 

T 

c 

0 

<1 

CVJ 

(0)/kTc 

per  atom 

The or . 

Exptl 

.  The or 

.  Exptl. 

Pb 

4 

1.38 

7.20 

7.22 

4.45 

4.43 

Pb0.80T10.20 

3.8 

1.23 

6.96 

7.08 

4.10 

4.04 

Pb0.60T10o40 

3»6 

1.02 

6.22 

6  „  16 

3.78 

3.81 

Pb0.40T10.60 

3  o4 

0 . 68 

4.51 

4.60 

3.51 

3.43 

Pb 

4 

1o38 

7.20 

7  0  22 

4.45 

4.43 

Pb0 .75InO .25 

3.8 

1.29 

6.92 

6.86 

4.31 

4.36 

Pb0.60In0.40 

3.6 

1.21 

6.68 

6.63 

4.19 

4.22 

Pb0.40In0.60 

3o4 

1.18 

6.65 

6 . 60 

4.11 

4.16 

For  the 

Pb-Tl  case,  Wu 

explains  the  : 

reduction  in  the 

ratio  2AQ(0)/kTc  in  the  following  way.  As  the  electron 
concentration  decreases,  the  effective  ionic  charge  also 
decreases,  as  does  the  Fermi  momentum.  Hence  the  momentum 


’ 


\ 

• 

• 

«• 
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space  available  for  the  integral  which  determines  the 
electron-phonon  coupling  constant  a  (co)  (equation  II- 38) 
is  also  reduced. 

For  the  Pb-In  case,,  although  the  electron  concentra¬ 
tion  decreases  along  with  the  electron-phonon  coupling 
strength,,  the  transverse  phonon  peak  shifts  to  lower  fre¬ 
quencies  where  it  enters  the  integral  for  the  energy-gap 
with  greater  weight.  The  two  effects  tend  to  compensate 
each  other  somewhat  and  the  ratio  decreases  more  slowly. 

In  attempting  to  explain  the  pressure  dependence  of 
2A  (.0)/kT  for  pure  Pb  along  these  lines,  one  notes  that 
although  the  electron  concentration  per  unit  volume  is 
slightly  increased  with  pressure,  the  phonon  modes  are 
shifted  to  higher  energies.  The  latter  dependence  tends 
to  reduce  the  coupling-constant  through  its  frequency 
dependence,  and  tends  to  reduce  A  (0)  by  decreasing  the 
effective  low  energy  modes  entering  its  calculation.  It 
is  of  interest  to  note  that  although  the  electron  concen¬ 
tration  increases  under  pressure  (and  the  Fermi  surface 
expands),  the  actual  density  of  states  decreases  in  the 
region  of  the  Fermi  level  (Garfinkel  and  Mapother  (1961)), 
so  that  it  is  hard  to  estimate  the  concentration  effect 
on  the  integral  of  the  coupling  constant  a  (go)  .  This 
effect  is  discussed  in  more  detail  in  section  VI. 5. 
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VI . 4  Changes  in  the  Thermodynamic  Properties  with  Pressure 
Garfinkel  and  Mapother  (1961)  have  carried  out  magnetic 
field  measurements  on  bulk  Pb  subjected  to  hydrostatic 
pressure „  They  report  the  following  results: 

(l)  The  pressure  dependence  of  the  zero  temperature  criti¬ 
cal  magnetic  field  is., 


dinH 

o 

dP 


(9.72  +  0.26)  x  10  bar  1 


(2)  The  pressure  dependence  of  txhe  critical  field  in  the 

limit  as  T  ^  T  is 

c 


dH 

(— ) 
dP 


-  (9.32  +  0.30)  x  10  0  gauss/bar. 


c 


(3)  They  use  result  (l)  in  estimating  the  pressure 
dependence  of  the  normal  metal  electronic  specific  heat 
which  is  given  by 


1  .  AS  *  y  1  ,  .  „  dH 

limo  -  =  7  T  =  (-)  T  =  -(-)  lim0  Hc(-^)p 
T-O  K  V  V  4tt  T-O  K  dT 


VI-15 


where  V  is  the  molar  volume,  and  y  is  the  coefficient 
of  the  normal  electronic  specific  heat  per  mole.  This 
is  done  by  expanding 


H0(P,T)  =  Hq(P)  +  A1(P)T2  +  A2(P)T4  .  .  . 

4  /  \ 

neglecting  terms  in  T  ,  and  estimating  A-.(P)  from  the 


VI-16 


l  * 


107 


temperature  dependence  of  H  .  This  gives 

(l/7*)(d7*/dP)  =  (l/Ho)(dHo/dP)  +  (l/A1)(dA1/dP) .  VI-17 


They  obtain  the  result 


* 

d  Jkiy  _  ^  , 

-  =  -(8.20  +  1.54)  x  10  J  bar  . 

dP 

*  p  2  /  v  /  v 

Since  y  =  2/3ir  k  N(0)  ,  where  N(0)  is  the  renormalized 
density  of  states.,  they  obtain. 


dinN(O)  ^  _n 

-  =  -(8.20  +  1 o54)  x  10  bar  1 . 

dP 


This  information  along  with  the  present  measurements 
on  dinAQ/dP,  can  be  used  in  estimating  the  pressure  depen¬ 
dence  of  the  condensation  energy.  As  mentioned  in  section 
II. 10,  Wada  has  shown  that  in  general. 


H 

—  =  N(0)  I 
8tt 


(where  I  =  Ja^  for  the  BCS  theory) 


VI- 18 

Here  I  is  a  function  of  the  renormalization  factor  and  of 

p 

the  complex  gap.  Introducing  the  ratio  l/(-|A^),and  using 
it  as  an  approximate  measure  of  the  coupling  strength,  one 
finds  that  for  P  =  0  the  ratio  is  0.83.  The  pressure  depen¬ 


dence  of  this  ratio  is  given  by 


din(l/&0/dP  =  2(dinHo/dP  -  dinAQ/dP)  -  dinN(0)/dP. 

VI- 19 


a 
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Substituting  the  relevant  measured  results  gives 

din(l/(fA^))  =  (8.7  +  3 o4)  x  10"6  barT1 

The  Indicated  error  includes  the  errors  for  dInHQ/dP  and 
d^nN(0)/dP  quoted  by  G-arfinkel,  One  finds  therefore,  in 
agreement  with  the  result  on  the  gap-Tc  ratio,  that 
increasing  the  pressure,  changes  the  properties  of  Pb 
towards  those  of  a  BCS  superconductor. 

The  pressure  dependence  of  the  energy  gap  can  be  used 


in  elec- 


-  x -  -  x  ov  >• — >  x. 

in  estimating  the  pressure  dependence  of  the  jump 
tronic  specific  heat  at  the  zero  field  s-n  transition e  As 
mentioned  in  section  II. 10,  Wada  has  used  strong- coupling 
theory  in  deriving  an  expression  for  this  jump 0  At  present, 
we  do  not  have  sufficient  information  available  to  carry 
out  a  comparison  with  his  theory,  however,  the  present 
work  can  be  compared  to  strong- coupling  model  calculations, 

p 

through  the  form  of  the  temperature  dependence  of  (A  (T)/kTp). 

The  BCS  assumption  of  an  energy  independent  gap,  leads 
to  the  prediction  that  the  electronic  specific  heat  should 
exhibit  a  discontinuity  at  T  given  by  the  expression. 


AC 


el 


3  d(£o(T)/kTc) 


7T 


2^ 


VI-20 


dt 


where  the  limiting  slope  is  to  be  taken.  The  BCS  weak- coup¬ 
ling  limit  result  for  this  is  1.43.  Figure  18  shows  a  plot 


. 
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of  (A  (T)/kT  vs  t  for  Junction  1345EQ  The  limiting 

O  O 

slope  for  the  P  =  0  curve  results  in  a  value  of 
AC^/7^  =  2.22  +  0.05.  The  same  slope  has  been  measured 
for  four  other  junctions.  One  can  also  estimate  values 
of  A-^(a)  =  0)  for  t  =  O.98  (figure  3)  and  t  =  O.96  from 
the  model  calculations  of  Scalapino,  Wada,  and  Swihart 
(1965)  and  Swihart,  Scalapino,  and  Wada  (1965).  If  these 
values  are  plotted  on  the  same  graph,  and  the  bulk  transi¬ 
tion  temperature  is  used.,  the  theoretically  predicted  slope 
is  also  2.2.,  which  indicates  excellent  agreement  between 
theory  and  experiment 0  If,,  however.,  one  takes  the  calori- 
metrically  measured  jump  in  the  electronic  specific  heat 
(Shiftman  et  al  (1963)*  Decker  et  al  (1958)*  and  Neighbor 
et  al  (1967))  and  makes  use  of  the  low  temperature  deter¬ 
mination  of  7  (van  der  Hoeven  and  Keesom  (1965)  and  N.E. 
Phillips  (1967))*  the  predicted  ratio  using  equation  VI- 20 
would  be  2.65  +  0.06. 

It  is  felt  that  this  discrepancy  between  the  present 
work  and  the  calorimetrically  determined  value  is  more 
apparent  than  real.,  since  the  framework  surrounding  the 
development  of  equation  VI- 20  does  not  really  apply  for 
strong- coupling  Pb ,  Swihart  (1962)  interpreted  the  discrep¬ 
ancy  between  the  BCS  value  and  the  calorimetric  value  of 
2.65  as  being  due  to  an  increase  in  the  gap  as  one  moves 
away  from  the  gap  edge.  This  is  probably  the  case,  and. 


, 

- 
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moreover,  this  argument  can  be  used  In  comparing  the 
present  work  to  the  calorimetric  result  since  choice  of 
A(cd  =  0)  in  forming  the  plot  of  (A/kT  )  vs  t  is  an 
approximation  to  the  energy  dependence  of  the  gap.,  while 
the  calorimetric  measurement  naturally  embodies  the  cor¬ 
rect  energy  dependence  of  A(cd)  . 

Although  Eliashberg  (1963)  may  have  been  the  first 
to  realize  the  importance  of  the  strong-coupling  electron- 
phonon  interaction  on  this  result,  Wada  was  the  first  to 
account  for  it  in  a  proper  way0  This  will  undoubtedly 
be  borne  out  when  a  direct  comparison  to  his  theory  can 
be  made  e 

From  the  shift  in  the  slope  of  (A/kT  vs  t  with 
pressure  (figure  18),  it  is  tempting  to  speculate  that 
the  specific  heat  jump  should  decrease  as  the  pressure 
is  increased.  This  should  again  be  explainable  using 
Wada’s  theory  and  the  pressure-dependent  phonon  spectrum 


for  Pb . 


* 

. 

' 

c.  >r 

v 
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VI . 5  Other  Shifts  In  Pb  Under  Pressure 

In  equation  11-42 *  the  density  of  states  (conduc¬ 
tance)  was  expanded  as 

Nm(ffl)  A^(cd)  -  Ap(cil) 

a  =  — -  =  1  +  - p— 1 -  +  .  .  .  .  VI-21 

N(0)  2a> 

If  an  energy  region  is  considered  where  A^  ~  0,  then  one 
gets 


d^n(a  -  l) 


dP 


dinA^(a)) 


a)=const 


dP 


VI-22 


a*=const 


This  condition  is  satisfied  for  the  gap-edge  region  of  the 
density  of  states.  Pressure-induced  shifts  in  the  conduc¬ 
tance  at  constant  energy  in  the  region  near  the  gap  edge 

are  listed  in  Table  9  for  2°K,  4.2°K,  and  for  t  -  T/T 

{  c 

Q 

close  to  lo  Typical  2  K  shifts  are  shown  in  figures  20 
and  21,  while  figure  23  shows  the  pressure  amplitude 
shift  for  t  -  0.92  and  t  =  0.95. 

Structure  amplitude  reductions  were  noticed  in  the 
region  of  and  when  the  junctions  were  subjected  to 

a  pressure.  These  effects  can  be  seen  in  figures  14,  15, 
and  16.  The  structure  is  associated  with  the  phonon 
emission  process,  and  should  be  sensitive  to  a  reduction 
in  the  electron-phonon  coupling  strength.  Defining  the 
structure  amplitude  as  S  =  |  a  -  1  |  ,  for  both  the  trans¬ 
verse  and  longitudinal  regions  (see  figures  15  and  16) , 
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TABLE  9 


Temperature 

Junction 

1 — 1 

1 

i — 1 

1 

t> 

0 

6P  din(cr-l)/dP 
(bar)  ( 10^/bar) 

2°K 

13450 

-0.037/0.765 

3172  o0 

-1.52 

o 

CV1 

1345c1 

-0.040/0.765 

3447.5 

-1.52 

4.2°K 

1345B 

-0.040/0.700 

2551.2 

-2.24 

4.2°K 

1345s 

-0.020/0.700 

1309 0 9 

-2.18 

t  =  0.95 

1345  A 

-0.30/7.200 

2378.8 

-1.75 

t  =  0.92 

1345A 

-0.50/11.80 

2378.8 

-1.78 

t  =  0.95 

1246a 

-0.20/7.200 

2137.3 

-1.30 

t  =  0.92 

1246a 

-0.40/11.80 

2137.3 

Average 

-1.59 

=  -1.73 

one  is  again  tempted 

dependence  of  the  gap 

to  associate  this  with  the 

through  the  relation. 

pressure 

din  |  o 

-  1  1 

dinS 

dinAAv 

VI-23 

dP 

<n=const  dP  co=const 

dP 

a)=const 

An  average  value  of  A 

^  is  taken,  since 

there  is 

no  way  of 

p  P 

estimating  dinA^(cu)/dP  and  dinA~(<D)/dP  separately  in 


this  energy  region.  Several  values  of  pressure  induced 
shifts  in  S  at  the  longitudinal  and  transverse  phonon 
energies  are  listed  in  Table  10. 


«i 

TABLE  10 


JCN 

<5S(a>t) 

esnp 

6P 

dinS(co.  ) 

1/ 

dlnS(o) .) 

Xj 

S(cot) 

S(oj^) 

dP 

dP 

(bar) 

( 10 //bar) 

( 10^/bar ) 

1345c 

-0.60/10.00 

-0.70/10.00 

3447.5 

-1.72 

-2.03 

1345D 

-0.50/9.90 

-0.50/9.60 

2757.2 

-1.83 

1 

1 _ 1 

CO 

00 

1345c1 

-0.60/10.00 

-0o 60/9.80 

3172.0 

-1.89 

-1.93 

Average 

-  -1.8l 

-1.95 

Taking  an  average  of  the  structure  shifts  in  Tables  9  and 
10  gives  as  a  final  average, 

d  in  (a  -  l)  _R 

-  =  -(1.81  +  0.35)  x  10  D  bar  ,  VI-24 

dP 


Using  VI-23 9  this  gives. 


din  A 


dP 


Av  =  -(9.1  +  1.7)  x  10  6  bar  1 


Since  the  directly  measured  gap  shift  was 


din  A  r 

- -  =  -(9.9  +  0.8)  X  10 

dP 

(see  section  VI.3)*  good  agreement  between  the  two  Indepen¬ 
dent  results  is  obtained. 

* 

The  pressure  dependence  of  y  as  obtained  by  Garfinkel 
and  Mapother  (1961)  provides  information  about  the  magni¬ 
tude  of  the  pressure  derivative  of  the  renormalization 
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parameter  Z  (0).  Writing  the  renormalized  density  of 
states  as  a  product  of  the  free  electron  expectation 
times  Z  (0)*  one  has,, 


* 

7 


2tt2  p 

-  kdN(0) 

3 


k2Nf.ei°)  Zh°)' 


VI-25 


where  N(0)  always  refers  to  the  density  of  states  per 
unit  volume.  Therefore,, 


diny  dinNf  (0)  dfnZ  (0) 

— . -  =  - LiHz -  +  - 2 -  .  VI-26 

dP  dP  dP 


The  free  electron  model  predicts  a  pressure  dependence 

dfnN.  (0) 

- -  -  +  k/3  VI-27 

dP 

where  K  is  the  compressibility  and  the  density  of  states 
per  unit  volume  is  used.  Anderson  and  Gold  (1965)  and 
Anderson.  0l  Sullivan.,  and  Schirber  (1967)  have  found 
that  the  Fermi  surface  for  Pb  scales  with  pressure  at  a 
rate  close  to  the  free  electron  prediction.  Therefore., 
using  K/3  =  0o68  x  10~6  bar"1  for  dinKL  (0)/dP,  and 
d^ny  /dP  =  -(8.20  +  1.53)  x  10  ^  bar-1.,  due  to  Garfinkel 
and  Mapother  (1961)*  gives 


dfnZ  (0)  r  .. 

- - -  =  -(8088  +  I.53)  x  10  b  bar  .  VI-28 

dP 
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From  McMillan  and  Rowell  (1965)  we  have. 


Z  (0) 
n\  / 


m 


00 

1  +  2  / 
o 


(co)F(o))do) 


=  2.33  +  0.02. 


VI-29 


Therefore, 

dinZn(0)  dinm* 
dP  dP 

and  the  effective  mass  is  decreased  toward  the  free- 
electron  value.  This  is  again  in  keeping  with  the  idea 
of  the  strong-coupling  superconductor  becoming  more  like 
a  weak-coupling  BCS  metal. 

This  information  can  be  used  in  estimating  the  mag¬ 
nitude  of  the  pressure  derivative  of  the  electron-phonon 
coupling  constant  a  (cd)  .  If  one  assumes  that  da  (cD)/dP 
is  independent  of  energy,  then  from  VI-29* 

•X* 

dm  p  00  F(cD)doD  00  F(co)dcD 

-  =  2a“d(  /  -  )  +  4a da  / - .  VI-30 

m  O  CD  O  CD 

Now, 

00  F(aD)dcD  9Nco2dcD 

din(  /  -  )/dP  =  din(  /  - )/dP  VI-31 

O  CD  %  O  CD  CD^ 

-dincD^ 


dP 


VI-32 


■ 


vnllQt'oo 


which  can  be  represented  as 


din(  / 
o 


F(o3)  do3 


03 


)/dP 


-din0D 

dP 


where  3N  is  the  number  of  degrees  of  freedom,  and 
is  chosen  to  be  a  Debye  phonon  spectrum.  We  have, 
therefore. 


dm 


m 


* 


,2 


=  -2a  (din0D)  / 


00  F  (  03  )  do3 


+  4ada  / 


00  F  (  03  )  d03 


O  03 


O  03 


00  F(03)d03 

I - 

O  03 


[4ada  - 


2a  din0D] 


* 


ml  p 

( - 1)  — p  [4ada  -  2a^din0n] 

m  2a^ 


Multiplying  by  m/m  gives 


dinm  =  (l  -  m/m  )[2dina  -  din©^]  . 

For  m/m  =  1/2.33 

*  2 
dinm  dina  din©n 

-  =  0.571  [ - -  ]  . 

dP  dP  dP 

Taking  dinO^/dP  as  an  average  between  dino9_/dP  and 
dino3^/dP  results  in 

din©/dP  =  (6.0  +  1.3)  x  10~6  bar'1. 
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VI-35 

VI-36 

VI-37 

VI-38 


VI-39 


(  -  - 
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Using  the  value  of  dfnZn(0)/dP  =  d^nm  /dP  from  equation 
VI-28,  results  in 

dina2  ^ 

-  =  -(9.6  +  2.2)  x  10  bar  .  VI-40 

dP 

Thus  the  electron-phonon  coupling  constant  decreases 
with  pressure  at  the  same  rate  as  the  effective  mass 
(renormalization  factor)  does.  Interestingly,  this  is 
approximately  equal  to  dfnAQ/dP  and  one-half  of 
d in (cr  -  l)/dP.  Therefore  one  can  tentatively  assign  the 
pressure  dependence  of  the  conductance  structure  with  a 
squared  dependence  in  the  electron-phonon  coupling  con¬ 
stant.  That  is. 


di!n(a  -  1)  2dinA/w  2dina2 
dP  dP  dP 


VI-41 


or, 

(a  -  1)  oc  A2v  cc  (a2 ) ^ .  VI-42 

2 

Although  the  estimate  for  dina  /dP  is  admittedly 
crude,  it  seems  hard  to  believe  that  the  interaction 
strength  increases  with  an  increase  in  pressure,  as  has 
been  suggested  by  Garfinkel  and  Mapother  (1961).  Using 
the  BCS  expression  for  the  condensation  energy  (a  poor 
assumption  for  strong-coupling  Pb),  they  predict  that  the 
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interaction  strength  V  should  increase  with  pressure. 

They  then  conclude  that  Tc  decreases  with  pressure 

because  the  change  in  the  density  of  states  over-rides 

the  change  in  V.  Although  a  rigorous  comparison  between 
o 

dinct  /dP  and  dinV/dP  cannot  be  carried  out  —  and  indeed 
is  not  justified  in  this  case  because  of  the  crudeness 
of  both  interpretations  —  it  is  intuitively  more  rea¬ 
sonable  to  assume  that  the  reduction  in  T  (and  also  in 
the  gap-Tc  ratio)  is  associated  with  the  interaction 
tending  toward  a  BCS  weak-coupling  limit. 

An  article  by  McMillan  (1968),  which  applies  to  the 
results  in  this  section,  has  recently  appeared.  Accord¬ 
ingly,,  this  discussion  is  continued  in  Appendix  B. 
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VI . 6  The  Recombination  Effect 

The  recombination  effect  as  proposed  by  Scalapino, 
Swihart,  and  Wada  (1965)*  was  outlined  in  section  II. 9. 

The  effect  arose  theoretically  in  the  finite  temperature 
calculations  of  A-^(cd)  and  A^(cd)  for  Pb.  In  particular, 
for  temperatures  close  to  T  ,  a  new  negative  peak  in  A^cd) 
at  energy  go  +  A  (T)  =  gu-^  was  founds  'where  gd^  is  the 
energy  of  a  predominant  phonon  group.  This  new  peak  in 
A~(gd)  was  reflected  as  additional  structure  in  the  tunnel 
density  of  states,  which  they  then  calculated  for  t  =  0.92 
and  t  =  0  095 .  The  process  ostensibly  involves  combination 
of  an  injected  electron  with  a  thermally  excited  quasi¬ 
particle  such  that  a  phonon  is  emitted  and  the  pair  then 
becomes  a  ground  state  pair.  The  process  Is  most  probable 
when  the  phonon  energy  equals  the  predominant  transverse 
mode  energy  ao^ . 

Khanna  and  Woods  (1966)  have  investigated  this  effect 
in  Pb  with  only  limited  success.  They  find  experimental 
tunnel  conductance  curves  which  should  correspond  to  lower 
reduced  temperatures  from  the  theoretical  picture,  and, 
at  the  same  time,  their  results  do  not  seem  to  scale  with 
temperature.  Thus,  further  studies  of  the  effect  were 
undertaken  as  part  of  this  work. 

Since  the  conductance  amplitude  is  highly  temperature 
sensitive  in  this  temperature  region,  accurate  temperature 
control  and  calibration  are  required.  The  thermometer 
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calibration  accuracy  is  estimated  at  0.1$,  while  the 
conductance  circuitry  is  accurate  to  0.2$.  This  results 
in  an  estimate  of  the  uncertainty  in  the  tunnel  conduc¬ 
tance  of  approximately  +0.003  for  the  gap-edge  region. 

Figure  22  shows  the  tunnel  conductance  (density  of  states) 
for  four  reduced  temperatures  near  t  =  1  obtained  from 
junctions  1345  and  1345A.  Shown  as  well,  is  the  conductance 
predicted  by  Scalapino  et  al  for  t  =  0.92  and  t  =  0.95. 

The  average  magnitude  of  the  conductance  is  in  good  agree¬ 
ment  with  the  predictions.  However,  the  fact  that  the 
experimental  conductance  exceeds  the  calculated  value 
between  approximately  1.7  and  3.9  meV  is  direct  evidence 
that  the  peak  in  A^cd)  is  smaller  than  calculated.  A 
reduction  of  Ap(a>)  in  equation  VI-21  would  increase  a  at 
the  gap  edge  and  would  also  result  in  a  reduction  of  the 
predicted  structure  at  V  =  cd  =  -  A  (T)  .  The  zero- 

temperature  tunneling  results  on  Pb  by  McMillan  and 
Rowell  (1965)  shows  that  the  transverse  peak  is  broader 
than  assumed  by  Scalapino  et  al  in  their  model  calculation, 
so  that  use  of  a  more  realistic  lattice  spectrum  would 
lead  to  a  more  shallow  minimum  in  A^co)  and  consequently, 
better  agreement  with  the  experimental  results. 

The  interpretation  is  further  complicated  by  the  fact 
that  the  phonon  emission  anomaly  near  cd^_  +  A  (T)  is  not 
sufficiently  separated  from  the  recombination  anomaly 
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Figure  22 


Normalized  conductance  a  for  junctions  1345 
and  13^5A  (full  lines);  the  parameter  gives 
the  reduced  temperature 0  The  dashed  lines 
are  calculated  normalized  conductance  for 
t  =  0.92  and  0o95o  The  arrows  indicate  the 
energies  -  A  (T)  and  co^  +  Aq(2). 
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near  a>,  -  A  (T)  for  t  close  to  1  where  A  (T)  is  small. 

It  is  possible.,  therefore,  that  the  recombination  effect 
is  entirely  absent,  Ir.  order  to  test  for  it  experimen¬ 
tally,  some  process  with  better  energy  resolution  than 
that  of  tunneling  at  this  temperature,  would  have  to  be 
used , 

Figure  23  shows  a  shift  in  the  conductance  amplitude 
under  pressure.  This  is  thought  to  be  associated  with 
changes  in  the  electron-nhonon  coupling  constant  as  dis¬ 


cussed  in  section  VI 0 5 


. 


Figure  23 


Normalized,  conductance  a  for  junction 
1345A  at  reduced  temperature  t  =  0.92 
and  t  =  0.95;  and  for  pressures  P  =  0  b 
and  P  =  2393.0  barc 
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CHAPTER  VII 

CONCLUSION 

The  normal  tunnel  resistance  of  Pb-Al  junctions  was 
found  to  be  reversibly  dependent  on  pressure.  The  pres¬ 
sure  dependence  Is  associated  mainly  with  compression  of 
the  oxide  layer  and  results  roughly  in  an  exponential 
dependence.  This  suggests  the  use  of  tunnel  junctions 
as  four  terminal  pressure  monitoring  devices. 

The  pressure  dependence  of  the  transition  tempera¬ 
ture  of  the  Pb  film  was  found  to  be 

dT  /dP  -  -(3.60  +  0.35)  x  10~6  °K/bar, 

0 

in  agreement  with  values  obtained  on  bulk  Pb .  The 
Gruneisen  y’s  for  the  phonon  peaks  in  the  high  dispersive 
region  of  the  phonon  spectrum  in  Pb  were  measured  at 
7t  =  2.52  +  0.74  and  =  3.32  +  0.49. 

The  zero-pressure  recombination  effect  as  proposed 
by  Scalapino,  Wada,  and  Swihart  (1965)  was  tested.  Re¬ 
sults  indicate  that  the  effect  may  be  due  to  the  rather 
idealized  choice  of  the  Pb  spectrum  used  in  the  calcula¬ 
tion.  The  temperature  variation  of  the  energy-gap  as 
obtained  from  zero  conductance  data  was  found  to  follow 
the  BCS  relation  to  within  a  few  percent,  in  agreement 


' 
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with  the  theoretical  predictions  of  Scalapino  et  al. 

The  shift  in  the  gap  parameter  Aq  with  hydrostatic 
pressure  was  measured  for  the  first  time.  The  average 
value  obtained  was 

dinAQ/dP  =  -(9.9  +  1.0)  x  10  ^  bar  2 . 

This  leads  to  a  pressure  dependence  in  2AQ/kTc  given  by. 


din(2A/kT  )  N  _6  .1 

- 2 - 2_  =  -(5.3  +  0.9)  x  10  bar  1 . 

dP 

t=const . 

The  pressure  dependence  is  thought  to  be  associated  with 
a  change  in  the  electron-phonon  coupling  constant,  such 
that  the  superconductor  tends  toward  a  weak- coupling  BCS 
behaviour  under  pressure c  Further  evidence  to  this  effect 
is  provided  by  the  fact  that,  under  pressure,  the  conden¬ 
sation  energy  tends  toward  a  BCS  prediction. 

The  conclusion  by  McMillan  (1968)  that  the  parameter 
A  defined  by 

0 

00  CL  (co)F(d))da) 

A  =  2/  -  , 

O  CO 


is  dependent  on  frequency  in  the  manner. 


A 


1 

<co2> 


has  been  verified.  Within  this  approximation  and  using 


■ 


f  . 


%  , 
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dinZ(O)  d in ( 1  +  A) 

dP  dP 


This 


dinZ(O) 


-(6. 


70 


+  1 


35) 


x  10 


-6 


bar 


-1 


dP 

is  in  good  agreement  with  the  value 
dinZ(O)  din7  dinN«  (0) 

 X  #  # 

dP  dP  dP 

-  -(8.88  +  1.53)  x  10" °  bar'1. 


obtained  using 

dirry*  c  _n 

-  =  -(8.20  +  1.53)  x  10  bar  1 

dP 

due  to  Garfinkel  and  Mapother  (1961),  and  dinNf  „  (0)/dP 
per  unit  volume  given  by  K/3. 
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APPENDIX  A 


Pressure  Shifts  In  F(gd) 

Consider  first  an  energy- independent  Gruneisen  y 

2 

and  a  constant  electron-phonon  coupling  constant  cl9  in 
the  region  of  an  and  go  0 


dfnco 

dinV 


A-l 


Under  pressure  P  each  frequency  changes  such  that 
0)f  =  GO^  +  Aon 

where  go^,  and  an  are  final  and  initial  frequencies. 

Writing  -  Aon,  and  remembering  that  the  compressibility 
is  given  by 

K  =  -  _ 

V 

we  have,, 


dV 

(— )T 


Ago/go  P 
AV/V  P 


A-2 


so  that., 

Ago  =  KPyoo  A- 3 
and.,  A  =  1  +  KPy  A-4 
Writing  the  shifted  frequency  spectrum  as 


F(a>)  =  AF0(o>/A) 


A-5 


and  remembering  that  the  total  number  of  degrees  of 
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freedom  is  3N, 

00  00 

/F  (ai)doi  =  3N  =  jF(oj)da). 
o  o 

We  have  also 

00  CO 

J  AF  ( cd/A)  das  =  /  FQ(a))dcD, 
o  o 

and 

00  00 

/  AAFq  ( a^/A )  dcp  =  J  Fq  ( a))  dcD^ 
o  A  o 

so  k  =  1  3  whence  F(cd)  =  IF  (co/A)  . 

A  A  0 


A-6 


A-7 


A-8 


A-9 


Differentiating., 


dF  (  cjd) 
den 


1  dF  ( a)  / A ) 
A2  den 


A- 10 


If  03  is  the  location  of  a  maximum  in  F  (gd)  such 
m  ox  ' 

that  F^(gd)  =  0}  then  os  =  Acom  is  the  location  of  a 
maximum  in  F(gd)  . 

Therefore  the  shift  in  the  maximum  gives  A  from 
which  one  gets 

7  =  A  -  1  .  A-ll 

KP 


We  now  consider  the  case  of  an  energy-dependent 
7.  If  we  assume  the  coupling  constant  is  slowly  and 
smoothly  varying  with  pressure  near  and  but 

is  essentially  independent  of  energy  in  either  the 
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transverse  or  longitudinal  branches,,  we  introduce  a 
pressure  dependent  term  so  that 

F(co)  =  A (  P)  F  (oj/X(co)  )  .  A- 12 

Differentiating  with  respect  to  energy, 

F'(o))  =  A  (  P)  F^oj/Afa))  )  )  .  A-13 

A  A2  da) 

Again  for  a  maximum  in  F  (a>/A(a)))  such  that  Fq(o)/A(o>)  ) 
equals  zero,  a  maximum  will  occur  in  F(cd)  such  that 

a)  =  A(cjq)cjq_,  and  y(a))=A(a))-l.  A-l4 

m  KP 

Note,  however,  that  when  A  (as)  =  03A'(cd)  that  F’(o))  =0 
as  well.  These  new  extrema  in  F(oj)  are  associated 
with  extrema  in  dA(a>)/da).  These  will  probably  be  of 
the  same  relative  importance  as  the  energy  dependence 
of  shifts  in  the  electron-phonon  coupling  constant  with 
pressure  which  we  ignored  in  our  original  assumptions. 

It  is  safe  to  assume  that  the  conductance  structure  shifts 
associated  with  these  secondary  effects,  will  be  much 
smaller  in  amplitude  than  the  structure  changes  in  the 
predominant  and  co^,  modes.  If  this  is  the  case,  then 
we  will  not  detect  them  anyway. 


* 
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APPENDIX  B 


A  Recent  Article  on  Strong- Coup ling  Theory 

Since  the  completion  of  this  thesis.,  a  publication 
by  W.L.  McMillan  (1968)  has  appeared.,  regarding.,  among 
other  things,  the  frequency  dependence  of  the  electron- 
phonon  coupling  constant.  In  the  paper  McMillan  uses 
existing  strong-coupling  theory  along  with  experimental 
results  to  consolidate  the  known  information  about  both 
elemental  and  alloy  type  superconductors. 

He  defines  a  dimensionless  "electron-phonon  coupling 
constant"  as 


CD 

o 


A  =  2/  a 


2 


o 


(CD  )  P  ( CD  ) 

v  q  •  v  q ' 


B-l 


and  from  equation  VI-24,  it  can  be  seen  that 


* 

m 

zn(°)  =  —  =  1  +  A, 
m 


B-2 


so  that  A  =  1.33  for  Pb . 

Although  most  of  the  paper  is  concerned  with  bcc 
transition  metals,  some  of  his  findings  are  applicable 
to  Pb „  In  terms  of  A  defined  above,  he  obtains  an 
expression  for  the  electronic  heat  capacity  coefficient, 
7,  given  by 


•* 

7 


7 


27T 


2 

k^N 


bs 


(0)  (1  +  A). 


B-3 


v 


3 


'•  ' 
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Here  N,  (0)  is  the  "band  structure  density  of  states" 
per  unit  volume  (containing  the  Coulomb  interaction 
between  electrons)  and  1  +  A  is  an  enhancement  factor 
due  to  the  electron  phonon  interaction.  Clearly.,  the 
band  structure  density  of  states  in  B-3  Is  to  replace  the 
free  electron  density  of  states  in  VI-35  and  1  +  A 
takes  the  place  of  Z  (0) .  Provided  the  band  structure 
density  of  states  has  approximately  the  same  pressure 
dependence  as  the  free  electron  density  of  states.,  equa¬ 
tion  VI-28  should  still  hold. 

McMillan  finds  that  the  frequency  dependence  of  A 
for  metals  such  as  Al,  In,  and  Pb  is  empirically  given 

by* 

Constant  K 

A  ^  ~  — p-  B-4 

M<co  >  <co  > 

where  M  is  the  mass  and  K  is  a  constant.  Using  B-2  and 
B-4  one  has, 

*  K 

dinZn(0)  dinm  dln(l  +  A)  djn(l  +  2^ 

dP  dP  dP  dP  B-5 

where  (from  equation  VI-28), 

* 

dinm  _r 

-  =  -(8.88  +  1.53)  x  10"  bar  . 

dP 

The  pressure  shifts  of  co^  and  aj^  can  be  used  in  estimating 


r... 
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diki(l 


_K _ 

<a)2>^ 


dP 


From  Tables  3  and  4  one  has  the  following  values 


a)  (P=0)  =  4.585  meV 


CD 


g(P=0)  =  8.865  meV 


and  from  the  average  values  of  6cd^/6P, 


=  0.075 


for  SP  =  3172.0  bar 


6cd^  =  O.I63 


for  6P  =  2757-2  bar. 


Since  1  +  A  =  1  +  k/<a)  >  ^  2.33  (McMillan  and  Rowell 

(1965)) 


6(1  +  K/ojp/6P  = 


2.2873  -  2.3300 
3172.0 


bar 


-1 


K 


K 


din(!  +  <-2>)  1 

o 

dP  2.33 


^  +  <^2>) 
6P 


-5.78  x  10  6  bar  2 


Similarly 


6(1  +  K2>p/6P  = 


2.281  -  2.330 
2757-2 


ba] 


-1 


and 


K, 


din(l  +  —2) 
_ 

dP 


-7.63  x  10  ^  bar  2 


B-6 


Taking  an  average  of  the  two  gives 


. 
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din(l  +  A)  &  ^  05  *  _6  _1 

-  =  - —  -  -(6.70  +  1.35)  x  10  bar  \ 

dP  dP 

B-7 

where  the  error  is  estimated  at  about  20^. 

The  agreement  with  the  independently  measured  value 
of  dinm  /dP  by  Garfinkel  (VI-28)  is  good  considering  the 
assumptions  made,  and  thus  the  prediction  of  an  inverse 
squared  frequency  dependence  for  A  is  borne  out.  One 
can  also  relate  the  pressure  dependence  of  A  to  the  pres- 
sure  dependence  of  a  (co)  as  shown  below.  By  definition 

^o  a^  (co)F(cjo)  dco 

A  =  2  / - 

O  0) 

where  05  is  the  phonon  end  point  energy. 

2 

If  we  assume  a  (a>)  is  independent  of  energy,  then 

dA  d  “o  a2(co)F(ca)dco 

—  =  2 — (  /  - )  „  B-8 

dP  dP  o  05 

Using  the  development  of  equations  VI-30  to  VI-36  and  the 

•)£ 

assumption  of  a  Debye  phonon  spectrum,  A  replaces  m  /m  -  1 
and  one  has 


•X*  Q 

dA  m  dtox  din9n 

—  =  ( - l)[ - -].  B-9 

dP  m  dP  dP 


Again  assuming  the  renormalization  parameter  can  be 
written  as 

* 

m 

zn(°)  =  —  =  1  +  A, 

m 


. 


* 


m 

A  = - 1 

m 
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B-10 


and 


dinA  din«  din© 

-  =  [ - 


D 


]  . 


dP 


dP 


dP 


B-ll 


Substituting  from  VI-39  and  VI-40  gives 


din  A 


dP 


=  [-(9.6  +  2.2)  -  (6.0  +  1.3)]  x  10  °  bar 


-1 


-  -(15.6  +  3.5)  x  10  6  bar  -1. 


B-12 


But  from  B-4, 


din A  2 din© 


<\j 

— 


D 


dP 


dP 


~  -(12.0  +  2.3)  X  10  °  bar  1. 


B-13 


Thus  we  again  have  consistency  between  the  various 


results . 
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